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ABSIRACT 


The application of LORAN-C in the hyperbolic mod=2 as a 
positioning system for hydregraphic surveys was investi- 
gated. Observed LORAN-C time differences from a field test 
conducted in Monterey Bay, California were compared to 
calculated time differences determined from geographic posi- 
tions based on a microwave positioning system. Four methods 
were used to determine the calculated time differences. The 
first three methods weze (1) applying only the seawater 
Secondary Factor, (2) computing the time difference based on 
a Semi-Empirical TD Grid, and (3) applying ASF Correctors 
meom the DMAHTC LORAN-C Correction Table. The final me*xhod 
applied multiple observed ASF Correctors at five minute 
latitude and longitude intervals. By applying multiple 
Observed ASF Correctors, which was the most accurate method, 
eeoce Ss Meter 1 drms with a tane offset of 3 to 12 merters 
using the 9940 X-Y LORAN-C combination was obtained. Rased 
upon the results presented, it may be pessible *o uss 


LORAN-C for hydrographic surveys at scales of 1:80,000. 
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Lee RO DUCT LON 


mee USE OF LORAN-C 


In recent years there has been an increasing usage of 2 
LORAN-C receiver and the LORAN-C network as the primary 
horizontal control for such scientific studies as deep ocean 
dumpsites, marine fisheries studies [ Rulon, 1979], bathyme- 
tric surveys, and recently, @ EPeconpatssance hydrographic 
survey. Examples of bathymetric surveys conducted by the 
National Oceanic and Atmospheric Administration (NOAA), 
which have used Loran-C for positioning are: 

1) Su-100-7-79 Gulf of Alaska [NOAA H-9822, 1979], 

mio U-100-2-79 Guilt of Alaska [NOAA H-9823, 1979], 

ey >-D9OZ2Z—-WH=-82 U.S. West Coase | NOAA Ship Surveyor, 

6-2 6 

Recently an attempt was made tc use LORAN-C as’ the 
Pemnaing position control for a reconnaissance hydrographic 
Survey S-K902-WH-82 [ NOAA, 1982]. This was a special survey 
Bemaacted by the NOAA Ship Whiting in May 1982. Spececad 
Surveys are field examinations of very limited extent or 
scope and frequently require unique survey or atap Coli =¢— 
meer procedures (Umbach, 1976}. The purpose of this preject 
Meemeco Verity the existence and extent of reported shoaling 
Mee nace Sarecy fairways in the Gulf of Mexico [NOAA, 1982}. 

mne use of Loran-C aS a@ positioning system for basic 
hydrographic surveys has been very limited due to *«ne abso- 
lute accuracy cf the long range systsn. A basic hydro- 
graphic survey is defined as a survey which is so ccmplete 
that it need not be supplemented by other surveys. “It must 


be adequate *o supersede for charting purpeses all prior 





surveys" [Umbach, 1976}. Variables which affect the accu- 
racy of LORAN-C are signal propagation variations, weather, 
and sky waves. The effects of weather and sky waves on 
LORAN-C propagation are best described by Samaddar [ 1980] 
and the American Practical Navigator, {({DMA, 1977] respec- 
eave] Y. 

Signal propagation variations are due +9 the phase 
retardation of the signal as it passes over an all Seawater 


path, cver land paths, or partial seawater-land paths as 


compared to free space. Table I= summarizes phase retarda- 
tion changes [Mortimer, 1978). Errors due *~)> an all sea 
TABLE. £ 


Phase Retardation or Lag of Radio Waves 











— me so ae a eet 
Difference in Phase | 
i Representation Lag at 500 km | 
Prcpa da =2 on Compared with | 
Veloce s Wave in Line | 
i Propagation Path (kn /sec Abcve2 (a) | 
a 
| Vacuum 2997902. 5 = { 
| 
| 
| Direct wave through 299GI) | 170 
| earth's atmosphere 
| 
| 
| Ground wave over 299 560 220 
| S¢a Water | 
| 
| Ground wave over 293 692 Pou | 
{ rugged mountains 
ee ee ner Ie rm Sr cece cre eer ee ee eee a ee eed 


Meret Path are Known es che Secondary Factcr (SF) and erzers 
due +o a land path or mixed vath are known as the Additional 


memomaary Factor (ASF) (Speight, 1982}. 


10 





ASF Corrections in the LORAN-C system can be as large as 
plus or minus four microseconds, which is 600 meters on the 
baseline. In other areas with the same LORAN-C coverage, 
these corrections may be much larger due to the expansion of 
the dis=ance between adjacent hyperbolic lines of position. 
For example, at 329 N and 809 W, using lattice pair 9960-X, 
a four microsecond (sec) error will offset the 9960-X line 
of position approximately 2438 meters [Speight, 1982]. 


mee APPLICATION OF ASF CORRECTORS FOR NAVIGATION 


To compensate for the Loran-C positional errors caused 
Pee the ASF Correctors, the Secretary of ransportation 
tasked the Defense Mapping Agencv Hydrographic/Topogzraphic 
Center (DMAHTC), the National Oceanic and Atmospheric 
Administration (NOAA), and the United States Coast Guard 
(USCG) wien eweTOb Df “d=termining and applying «the 
Mga ioral Secondary Correctors for each Loran-C chain. 
Mas task was published in the Departmen* of Transportation 
Moty Nacional Plan for Navigation in the July 19, 1974 
Federal Register. Wess COErtect=ons should orovid= 95% 
assurance that a vessel could fix its pos 
NM) MWe thir che U. S. 
Coastal Confluence Zone (CCZ) and *h2 Great Lakes. Piemee 


Lz. One com ra PL=a— 


Meeed acctiracy of 1/4 nautical mile  ( 


tN 


is defined as: 


"the seaward approaches zo 1 she inner Dourdary cE 
Mec 2S he narbor entrance and the outer boundary of 
Moech =S 50 nautical miles offshore or the 3dge of «he 
Continental Shelf (100 fathom contour) whichever ‘is 
mecater. ' 
The 1/4-NM accuracy requirement alse affects the 
maumencal chart. The National Ocean Survey (NOS), which 


pueeeshes charts for the CCZ, engaged in a program with the 
MeeG end DAAHTC to provide the coastal navigator with charts 
Overprinted with lattices which mest 1/U-NM accuracy. The 


USCG, as operator ci the LORAN-C ctadionavigation systen, 


i 





conducts surveys to ensure that LORAN~-C coverage exists 
Peehon the CCZ and will be reponsible for the verification 
of 1/4-NM accuracy for all coastal LORAN-C service. Ea 
Bonjunction with NOS, -t assists in surveys cf coastal 
waters of the United States to allow production of LORAN-C 
charts based on observed field data to meet the standards 
set forth above [Speight, 1982]. 

Bmndee, foc LORAN=C civil need, prepares grid predic- 
rons from its data base. Based on analysis and verifica- 
tion of che predicted grid from a USCG and/or NOS survey, it 
Produces revisions to the initial grid predictions (Speight, 
io 2). At present, DUAHEC Nas provided NOS with ASF 
Corrected LORAN-C Lattices which are overprinted on the NOS 
Miaets., Each chart with ASF Correctors applied contains one 


of the following notes: 


“The LORAN-C lines _ of position overprinted cn this_chart 
have been repared for use with groudwave signals and 
Mem ecole) COMpPetsd-2=d only EOL g-neoretical propaga- 
tion delays, which have not yest been verified “by 
observed ata. Wome ae S aro YEant2oned not to rely 
Peale ly on the lattices in inshore waters. Skywave 


mee rections are not provided". 
a 


Meta LORAN-C ines of positicn overprinted on this chart 
have been prepared for use with groundwave signals and 
are emncuse Cec Toots vaga- On Gelaye Ccompuszed Eron 
observe data (ees 5 abe Galle Oned Anode co rely 
entirely on the lattices in inshore waters. Skywave 
MmemnecclLOns are not provided" (Speight, 1982}. 

Presently, omer e HemmNOS charrs - of 7:80,000 to 
120,000 scale covering the east coast, Gulf coast, and 
Great Lakes show LORAN-C lattices that have been compensatisd 
BeeeAdditzional Secondary Factors. Pact OL wie La=sices ‘on 
these charts have been constructed from DMAHTC data tapes 
“Mat provide adjusted LORAN-C readings for each rate at 
every five minutes of latitude and longitude. A tew 
lattices were constructed using a single ASF Correction for 


the entire chart area. Five minute data tapes were not 


qZ 





furnished by DMAHTC for constructing lines of positicn for 
LORAN-C rates on the West Coast Charts. OMe ciess Chares 2 
single average ASF Correction was used to adjust e¢éach 
lattice [NOAA, Marine Chart Division, 1982]. 

In addition to supplying corrected LORAN-C lattices for 
nauzical charts, DMAHTC prepares, distributes, and periodi- 
cally updates unclassified ASF LORAN-C Correction Tables 
Mmeperght, 1982]. MaewAse (GCE ECCtLI2 Lables are £Or preci- 
mon Navigation, utilizing digital computers to convert 
LORAN-C time differences to geographic coordinates [Speight, 
Moo 2}. Presently, the ASF correctors found in the LORAN-C 
Correction Tables were determined using theoretical propaga- 
tion delays. ASF corréectors listed in the tables are going 
to be updated with observed data and reprinted the first 
G@iacter of 1983 {Wallace, 1982}. 


Meee APPLICATION OF ASF TO HYDROGRAPHIC POSITIONING 


Schnebele [1979] investigated the possibility of using 
Loran-C as an electronic positioning system for hydregraphkic 
surveying. TomcOvuemldede. nae 10 “ont2tey Bay, California a 
Singie Additional Secondary Factor (ASF) appiied to offshore 
dines of position gave a root mean square error (drms) of 66 
meters for the West Coast 9940 Y-W pair and a predicted 42 
Meeert drms error for 9940 X-Y rates. 

The 42 meter prédicted drms is larger than Nelson's 
Sumeetee. Electric Co., 1979] findings in San Francisco Bay. 
He demonstrated, in adynamic mode, that the precision of 
LORAN-C was 60.8 meters 2 drms (30.4 meter 1 drms) with 2 
Memst Case of J/1i.2 meters 2 drms (35.6 meter 1 dims). A 
Mean difference or offset between the neasured time differ- 
ence and the calculated time diiference was 34 nanoseconds 


Meoreehe 9940-X rate and one nanosecond for the 9940-Y rate. 


fies 





He also obtained a precision of 38.0 meters 2 dims (19.0 
meters 1 drms) in the static mode. Nelson also states, that 
the above precision is only achievable if the user has a 
LORAN-C receiver which has the performance capabilities of 
those used in the experiment. The LORAN-C receiver must 
have "comparable Signal averaging time, extra notch filters, 


and attenuation of the signal" (General Electric Co., 1979]. 


tae OBJECTIVES 


The National Ocean Survey requires that hyperbolic 
control systems used for hydrographic surveying exhibit a 1 
drms of less than 0.5 millimetar at the scale of the survey 
Pmumbach, 1976 j. Although this requirement is generally for 
2mHz phase comparision systems, it can be inferred that i+ 


also applies to other hyperblic svstems such as LORAN-C. 


The scale routinely used for coast2 Surveys is between 
1340,000 and 1:80,000 (Umbach, 1976] yielding an allowable 
fmeeor Of 20 to 40 meters net includin SYS@=Na fen es rees. 


Beamebele [1979] concluded that hyperbolic LORAN-C, after 
meee vaing a2 Single ASF Corrector, =Seunsvieese ble gasor basic 
hydrographic surveying. 

Whether or not applying multipl2 Additional Secondary 
PactoOrs (ASS) to LORAN-C lines of position will zseduce the 
drms sufficiently to meet the accuracy standards set by the 


National Océan Survey Hydrographic Manual will be 


pv 
{in 


CS] 


U) 


Saaned in this study. Mhemcern WMuletiole ASF Cor=ector 
Memers tO the application of more than one corrsetor to 
LORAN-C lines of position over a given aréa. The variable 
BeeecOrrectors result from varying delays of the el¢ectromag- 
netic wave as it propagates over different land segments. 
Three methods cf applying multiple ASF Correctors were 
Hl 


tested. The first method was the application of a 
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Semi-Empirical Time Difference Grid Calibration Model devel- 
oped by The Analytic Science Corporation [ 1979]. The spon- 
soring agency was the United tates Coast Guard. The 
Semi-Empirical Model applies Secondary Paceors and 
Additional Secondary Factors for each geodetic position 
based on the distance over land, the distance over water, 
and the total distance using mean se2 water and land conduc- 
tivities. 

The second method which was investigated applies ASF 
Correctors found in the DMAHTC LORAN-C Correction Tables 
{DMAHTC, 1981] to LORAN-C lines of positions. These correc- 
tors were derived from the ground conductivities which have 
been determined in the field by a Toast Guard calibration 
team [U0.S. Naval Oceanographic Office, 1982]. 

Finally, a third method was pursued. ASF Correctors, 
which were determined by fleld observations, were applied to 
the LORAN-C lines of position. Thes2= were determined by 
computing the difference between “«h2 observed LORAN-C rates 
and the expected time difference which was calculated using 
Bout lines of positicn from a very accurate microwave posi- 
tioning system. These ASF Correctors were determined at 

£u 


meve Minute latitude and longitude intervals. 


is 





II. NATURE OF THE PROBLEM 


Meee CHE PRINCIPLES OF LORAN-C 


To understand the problems associated with LORAN-C when 


used during hydregraphic ®perations, One MUSeetzsse Undsr- 
Stand its principles of operation. LORAN-C is a low 
frequency, pulsed signal, hyperbolic, radio navigation 


system, employing time difference measurements of signals 
received by the navigator from at least three ground tran- 
smitting stations { Speight, 38 25 The stations ar? 
comprised of a master smanscient t ung station, two or more 
secondary transmitting stations which are strategically 
spaced several hundred miles apart and, if necessary, a 
Syeeem Area Moniter (SAM) Station {U.S. Coast Guard, 1974}. 

System Arsa Monitor (SAM) stations asscciated with each 
Menman-C chain apply differential-type correcticns +9 he 
rates in real-time. SAMMY Sceetons continuously monitor the 
Seamaes ©£rom all transmitters in the chain. T= the oberved 
Meee ricrence deviates by- more than 0.05 msec ‘from the 
expected value, then the appropriat? secondary adjus<s its 
emissicn delay time in order to remcve2 the error [Schnebele, 
79). 

The master and at least two secondary stations a4 
meeaeted Such that the signals from the <=ransmictting stations 
can be received =hroughout the desired coverage area. z 
Master station is designated by the letter "M" and the 
Secondary stations or slave stations are designatied W, X 
for 4 {.U.S. Coast Guard, 1974}. 

Mei Sta=i0nsS transmit on the common frequency of 100 
AZ interference between transmitters is avoided through 

2]. 


the use of time separation [Poppe, 19 After the master 
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station transmits a pulse, each secondary station delays its 
Own transmission for a fixed time, called the secondary 
coding delay. This ceding delay is synchronized through the 
use of cesium frequency standards at each station. The high 
stability and accuracy of these standards permit e¢ach 
station to derive it own time of transmission without 
Bererence to another stazion [DMA, 1977}. Secondary ceding 
delays are predetermined by system propagation times and 
equipment characteristics { Laurila, 1976}. 

The pulse from the master transmitter is distinguished 
from those of the secondaries through phase coding of the 
pulses. Phase coding refers t9 the inversion of «he nega- 
tive and positive peaks of the sine wave comprising the 100 
az Carrier portion of the pulse. The purpose of the phase 
Bed2ng is twofold: 


"FIIst, He pPermics automatic discrimination between the 
Master and _the various Seecodansy sta: -LONS, thereby 
ermitting all stations to be identified by their rela- 

ive timing with respect te the master' 
"Second, the phasec rend Deve d=sS EOteCeC tL On sagelns— 
excessively long skywave delays which would cause the 
late arrival ee the preceeding pulse to coincide with 
2G 


the pede ee edge oe ae nae DGnerton Of a pDUuLSe Db 
tracked" oppe, 


The signals are received by a mobile receiver where the 


differences in time of arrival of the master signal and 


Various secondary signals are measure and displayed cn the 
mem@rcator portion of the LORAN-C set. The accuracy of this 
time difference is increased by phase comparision “of the 


synchronized 100 kHz carrier within the master and secondary 
pulses" [Laurila, 1976}. This measured time difference (TD 
= in Microseconds) represents a hyperbolic Line of Posizion 
meee) [| U.S. Coast Guard, 1974]. The intersection of two or 
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more LORAN-C LOP's defines the 
femagure 2.1). When plotted on a 
the resultant hyperbolic MJlines 
Pestion { Speigh=, 1982}. 
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ber PHASE LAG 


In a vacuum, the velocity of radiated energy from an antenna 
for LORAN-C is 299792.458 km/sec. Since radiated energy 
cannot be shown pictorially, the ohase of the *ransmitted 
radiations is used. The lines of constan+ phase of the 
transmitted radiation are shown in Figure 2.2 by the curved 
lines labeled aa't, bb', and other similar designations. 
They define the wave front as it proceeds outward from thes 
antenna in all directions. The distance between each line 
of constant phase is one wavelenath (\) [Admiralty Manual, 
1965]. 


th 


H 
(D 


A (meters) = 299792.458 km/sec + 
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quency in kHz 
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Figure 2.2 Transmitted Radiation 


The velocity of the radiated energy in air de¢evends on 
d 


temperature, pressure, PUMLdeey,.. 2n 


ig 





surface over which the transmissions pass. The retardation 
of a transmitted wave is known as phase lag. When low 
frequencies are employed, such as LORAN-C at 100 kHz, the 
effects of change of temperature, pressure, and humidity are 
swamped by the effects caused by changes in the nacture of 
the surface over which the transmissions are traveling 
{Admiralty Manual, 19634). The General Electric Company, 
mem PO division, conducted a LORAN-C Signal Analysis 
Mepetiment under the directton of the U.S. Coast Guard. 
This experiment was conducted along the U.S. West Coast. 
The General Electric Company recorded an overall change of 
108 nsec and 116 nsec time of arrival from the master and 
X-secondary stations respectively after a storm [Samaddar, 
1980}. If ASF Correctors are as large as two microseconds 
Smeene West Coast [DMAHTC, 19817] then the ASF Corrections 
are 20 times larger in magnitude than weather effects for 
the 9940 West Coast LORAN-C chain. 


The change in transmission rat2s or phase lag are a 


(Dd 


result of the amount of e2nergy “«ransferred from ‘he “tran- 
emt -e-ead ~tadiation. This energy transfer depends cn the 
Mme Pticon dualities (inversely related to conductivity) of 
the surface over which they are passing and their wavelength 
fer frequency). The lower the conductivity and the longer 
the wavelength (or lower the frequency) the greater the 
transfer of energy, and vice versa. Seawater has a rela- 
tively high conductivity (5.0 mhos/meter). Land has a much 
Memwer conductivity, which varies from marsh (fairly high) to 
dry sand and rock {very low) [Admiralty Manual, 1965]. Me 
Seeeriant papers that discuss the electrical properties of 
Soil are those of Smith-Rose [1934] and Pressey, Ashwell, 
and Fowler [1956]. Smith-Rose [1934] found that the conduc- 


tivities for soil ranged from 0.18 mhos/mezer for 2 grey 


+ 


€ 


Clay with salt to .00001 mhos/meter for granits. 
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the phase lag increases with 
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Secondary Factor 


most intriguing Cee Was! a= 


interface where there is an extreme change in 


the 


land and sea. 


the 


ivites between Visualizing 


en three  ~dimensions, lower part of the 


Myethesa@teg of the ground, lags further and 


the upper part as the wave crosses the land. 
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coastline it suddenly encounters the much lower impedance of 
the sea, and in a very short distance the bottom of the wave 
tries to catch up with the top, as though the whole wave 
front were an elastic balloon. This is known as "phase 
recovery" [ Eaton, 1979]. Phase recovery was verified during 
tests on Decca ‘*ransmissions across the south coast of 
England by Pressey, Ashwell, and Fowler [1956]. 
mae deterxrMlGB@ation of the S€condary Factor for seawater 
Memerairiy direct since ‘the conductivity of Stawater (5.0 
mhos/meter) is fairly constant. But for land the conduc- 
tivity can vary depending on the type of soil and its water 
content [Smith-Rose, 1934}. 
Phase lag for radiated energy over land can be deter- 
Mined two ways: 
1) AsSign an average land conductivity +o the ASF 
Model. For example, the average conductivity for 
Ene SOil on *he west coast <s 0.003 mhos/meter. 


The average land conductivity will determine 


(D 


(iD 


eon 
average phase ratardation of the pane 


matey tsc SclLence Corporation ,1979]. 


2) Dereriienme GVeEry COnduccelyity for each portion 
of a line segment from the transmitter to the 
receiver. The mocal Or these conductivity 
segments constituting a land-water profile will 
determine the tctal phase retardation of the path 


fpspeight, 1982}. 


fee 2 De MODEL 


PoSitional fix accuracy Sing. GORAN=C wisS Ori marily 
= 


a 
dependent on a chart makers ability to accurately computs 


the expected difference in time-of-arrival (TOA) of received 


O 
G-Ourdwave Signals from the transmitting stations. Time 
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differences (TD), are the differences between the TOAS of 


the secondary and master transmitters. 


= TOA . - TOA (251) 
a x m 

= eee Staticn 

m= Master Station 


TOA computations are dependent upon an accurate knowledge of 
the signal vhase delay. 
The phase delay of a aqroundwave signal is generally 


expressed as: 
o = T + oe (252) 


= -- - SF 


where n is the surface refractive index, cC is the speed of 
light in a vacuum, R is th? range between che receiver posi- 
tion and the transmitting station. The primary phase delay, 
T, is the computed travel time of the LORAN-C pulse cover a 
Mestance equal to the transmitter-to-receiver great circle 


meemesength, taking into account the velocity of electromag- 


ct 


netic waves and the index of refraction cf the atmosphere. 
Mees secondary factor (SF) UemmerColrecc On =O “he pzimary 
phase delay and accounts for the phase laq. The dominant 
meen in (2. 2) is the primary phase delay (T). The SF is 
uSsuaily an order of magnitude smaller [The Analytic Science 
MemapOoration, 1979]. 
Thus, time-of-arrivals can be axpressed as: 
TOA, = T. ~ SF + cD. (233) 


TOA = z + SF (2.4) 
m mn mn 


}?- 


where CD is the true emission delay or coding jielay for th 
Meme —-C chain {The Analytic Science Corporation, 1979]. The 
fr 


coding delay is equel to a <zime delay plus a2 compuced on 
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way baseline time (Bc) which includjes the secondary phase 
correction for an ail seawater path. The oneway baseline 
time (Bc) is equal to the distance between the master and 
secondary transmitters in metérs divided by the propagation 
velocity of LORAN-C through the earth's atmosphere (299.691 
meters per microsecond [Navigation Department DMA, 1982)). 
Sea Table II for Ccding Delay values for th2 9940 chain 


TABLE ITI 


Coding Delay 9940 LORAN-C Chain 


| Pair 9940-w: CD + Bc = 11000 + 2796.90 = 13796.90 psec | 
9940-X: CD + Be = 27000 + 1094.49 = 28094.49 psec { 
PovO-vam cD + Be = 4OCO00 + 1967.2/ = 41967,27 psec | 








memmordan, 1979]. Combining equations 2.1, 2.3, and 2.4, ths 


(Dp 


Meage TD +S given by equation 2.5 [The Analytic Science 


Seeporation, 1979}. 


Picea T ) +) (Sr. = SF) + CD. (2.5) 
a as m = n ma 


fe SeUL-SMPLRICAL TD GRID CALIBRATION MODEL 


The Semi-Empirical TD Model was developed by The 
Analytic Science Capo bacon (1979 ] aa Reading, 
Massachusetts for the West Coast 9940 LORAN-C chain. 
sa@etar "time difference (TD) grid calibration techniques 
have been successfully employed to develop an accurate 
(approximately 100 nsec drms) calibrated grid for St. Marys 
Meyec LORAN-C chain", (The Analytic Science Corporation, 
09}. 
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1. Technical Approach 


The Semi-Empirical Model is based cn Millington's 
empirical approach for computing the secondary factor over a 
mixed (multiple-homeogensous segmen*) path which combines 
land and sea phase delays. The gensralized semi-empirical 
polynomial functional form for the SP of the LORAN-C station 
is given by: 

- 


ie 
eee = SF(T ,8.) = Maemo +) D costs). (2.76) 
j js jt J nee 3 


DAT. + 
oe ee ok 


where 
j = secondary (W,X, or Y) or master (M) station, 
nR- 
J 
Te ejciista t2On=- LO Use= paltaly phase delay, 
, e 
R = jth station-to-user gre2zt-circle path length, 
J 
8B = user path bearing angl= at the jth station, 
J 
K , K_ and tL aré vositive integers, 
1 2 
C.  andD are the station-dependent coefficients 
aL ape 
OL NarMonic tLe=ms in the nodel, 
- is the range-dependent coefficient of the model 
which may in general be station-dependent. 
Date Erom 27 coastal sites distributed along che 
? 
Meme COeSc and 122 land-sea sites distributed in the 


Souzhern California CCZ (between Point Arguesllo and San 
Diego - see Figure 2.4) were used in @ Kalman estination 
elgorithm to compute the uncertain coefficients of the land 
@paesca mModcls of the TD grid calibration algorithn. (An 
Mempranation Of Kalman filtering for *«he layman i 


S 
by Roger M. du Plessis [1967].) The vGal2zbreced 
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was used to compute TDs at eéch data site and the TD resi- 
duals (difference between measured and calibrated TDs) were 
examined. Adjustments were then made to the TD model struc- 
+ure in an attempt to further reduce the residuals. Phis 
process of adjusting the nodel structure is repeated until 
the residuals agree with the expected theoretical covariance 
associated with the TD model. The model which exhibited the 


"hest" performance was selected as *+he West Coast TD grid 


model. 
2. Generalized Range/Bearing Model 
The Generalized Range/Bearina (GRB) Model was 
selected as the "best" semi-empirical calibration mod¢l for 


the Wes* Coast chain. The semi-empirical function is: 


SF = AO + AIT + : (Gees niet DD cosis ) nsec (Ze 
j owls s i: 

where AO, Al, Cand D are the model coefficients, is the 
path beating angle measured positive clockwise from north at 
mermeme) (| W,X,Y or M) sStaticn and Tf is tha path range «so <=he 
men Station. The GRB model is relatively complex and is 
expected to exhibit superior performance. The 2xtzensive 
mcdel is based on knowing the distance overland (TL), he 


distance over water (TS), the total distance (T), and the 
Bae bearing angle 8 (Figtre 2.5). 

It was noted that the calibrated medel was expected 
memoe accurate ard applicable only over the extant of ranges 
end bearing angles embodied in «he calibration data. Hence, 
outside the region covered by the calibration data *he model 
May not be as accurate as within the data coverage region. 
Using the GRB model, a drms value of approexinately 0.8 pisee 
was expected in areas where land data alone was use ae 
calibrate the model. Inclusion of sea calibration data 
Begauced 2 dems value of 9.35 to 0.50 psec. 
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ELGure 2. 5 Mixed Dath [ID Geometry 


The time difference (TD) is expressed 
following equation: 


[eT ) t9(SP = SF) + CD, + Bb, psec 
a a mn a m i 
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where 


R e 
i: 
=o psec, 
a Cc 
nes 
m 
i) = — psec, 
m Cc 
R. = ith secondary station-to-user great-circle path length, 
a 
@. = ‘the SBSée of an 2lectromagnetic wave in a vacuum 
= 2.99792458 X 108 m/sec, 
nh = surfac2 refractive index 
= 1.000338, 
R = master station-to-user gr¢eat-circle path length, 


CD = ceding delay found in Table II, 
b = TD bias associated with the ith secondary 
station (psec) (Table III). 
Ag ps Ode 


MDs Ss 7 (b) -= isiec 





| TDW -0.854 | 
TDX ae 
| TDY Sess) 
u Oe eed 
m= On> [-s1 + S2 + 53 - Su + 55 + S6) (2.9) 


MWreeremle 5! 25 the SF Gof a Land vath of length T 


ae 

(nsec) from the jth station: : 

Beme = 0.795/T +¢+ 0.439 + (0.00245) T Sec. (23530) 
sj $7 

fhe terms SZ, $3, and S& combine =o make up <zhe 

Seeeomecrymneaceor for the Seawa=er path lengths. 
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TABLE V 


Land Coefficients 


i care 
| AO = 1.428 D 5 = 0.942 | 
x 
A1 = O00 156 c ; = Oa) 
C es 0.0 ee = 0.588 
wi y2 
= -0.711 D = 0.0 { 
( Ww y 1 | 
\ D = Osa D = 0.0 
| wi We | 
| D = 0.0 e = 1.010 | 
w2 m1 
| c = O20 EC = -0.196 
| x1 m2 | 
ec = 0.0 D = -0.893 
| X2 m 1 | 
| D = 0.0 D = = 05.3595 
Sl m2 \ 
ee ee a alae ete see tem ila Sin ce i ci icc Aewnse-mh es nh cmd 


ie DMAHTC MODEL 
1. Sea SF Model 


The equations for the Sea SF Model is: 


SF areeeweceee (Bo tT) usec, 2£ 1) S T s 537 psec, (2. 13a) 


ar 


met) + AZ #9 (A3 T) wsec, Lt T > 537 psec, (223 5) 


where T is the primary phase delay (or range) in microse- 
conds (usec); Ak and Bk (k = 1, 2, and 3) are the sea nodal 
coefficients used by DMAHTC in program MTDGRID [Funekoshi, 


1982]. The coefficients are found in Table VI. 
2- Land sf Model 


inemecoouu=s on used <5 resolye DMAHTC Tables ASF 


Meeeecticns is called Millington's Method fDMAHTC, 1981). 
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Sea SF Model Coefficients 


MCT ic eee eee Pomc oe 
deen Sedu oS 
A2 = -0.890758 

A3 = 0.00064576813 
iene ee Daradeee 
iso = =omonamoz 
|B3 = — 0.00032774815| 
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This methed is based on the premise that the phase distor- 
tion due to a composite land-sea path is the arithmetic 
average of the phase distcrtion found in the forward and 
reverse paths of the propagated sign2l (DMAHTC, 1981}. FOr 
example, in Figure 2.6 two azimuths have been drawn on the 
Map and are labeled as 2109 and 2359. Also placed on the 
map are the proper ground conductivities which have been 
determined in the field by the Coast Guard calibration <«+<¢an. 
meueecas Circle drawn on the appropriate chart or charts from 
mre LORAN-C Station ccordinates to the area under considera- 
tion spans various lengths of land and seawater. Baca 
length or segment will have a specific conductivity and 
Meeeance, The total of these conductivity segments, consti- 
m@e2ng a land-water profile, will determine the total phases 
retardation along that path [JU.S. Naval Oceanographic 
Meeace, 1982}. All azimuths and distances are compuzed 
based on the World Geodetic System (WGS) datum {[DMAHTC, 
eo 1]. The values of phas¢ fre 

Beweauctivity are tabulate a2 the National Bureau of 
Seanadards (NBS) Circular 573 (Speight, 1982}. 


Diemer OEM la mused SO derive the ASF Cotrection for 
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tascaee Ol LOL a daven ground 


the time difference for a master-siave transmitting s' 
Peatr is: 


ASF Correction = (-Slave Error) - (-Master Error) 
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Figure 2.6 Composite Land-Sea Path 
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presented in an azimuth array. Pnesudasdye:Ss a Series of 
geodetic azimuths radiating from the transmitter with 
correcticns computed at incremented distances along e¢ach 
azomuth. After the values are computed in the azimuthal 
array a cOmputer pregram rearranges them into a matrix form. 
The matrix form is the arrangement of corrections into rows 
and columns covering a specified geographical area at @ 
constant spacing. The purpose of the matrix is to enables 
+he corrections from two LORAN-C transmitters to be added 
algebraically, combined into a single matrix, and arranged 
in the desired tabular form. in=seraotlas Lorn 25 the body 
@eechne table (DMAHTC, 1981}. 


3. Table Description 


Fach table contains a complete chain. Figures 2.57; 
a0, aided ee depetct LORAN-C ASF Correctors for chains 
mero-W, xX, and yY for Monterey Bay, California. A table 
Secticn is prepared for each station pair (master station 
ee one slave station) in a LORAN-2 chain. As a rule the 
limits of the table coverage are determined by the range of 
the groundwave transmissions for the LORAN-C chain. Each 
page of corrections in the table covers an azea ‘threes 
degrees in latitude by one degree cf longitude, with correc- 
meeons printed in increments of five minutes of arc. Rats 


designation and page numbers are printed at th 


(p 


top cf each 
correction page. Those pages where latitude and longitude 
mete s Contain both land and sea are included but correc- 
Meees apply only fcr the area covered by the U.S. Coastal 
Menmciuence Zone (CCZ). Large land bodies and areas ousid= 
the CCZ are represented by blank spaces on the page. ASF 
correcticn values can be either positive or negative (pesi- 
tive values are shown without sign). AreaS requiring no 


correction show a zero value which in some cases in preceded 
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by a negative siqn indicating thet the zero results from <he 
reunding off of a value slightly less than zero (indicates 


mae trend or the correction). 
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F1gwme 2.8 9940-XK ASF Correction Tabis¢ 


Mie “fable can be entered directly by using the 
Ship's position determined to the nearest five minutes of 


arc in laczitude and longitude either by dead reckoning or 
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mequre 2.9 CGu0=¥ ASF Gorrection Tabl 


some other means. To find the page w 
menrect.on, the Page Indexes of the table should be 


utilized. These indexes show the limits and oage number of 
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all pages in the table. To locate the number of the page on 
which the desired correcticn is to be found zhe Page Index 
is entered with the ship's position. In some cases the 
ship's position will fall on the page limit in either lati- 
tude or longitude or both. These positions are repeated on 
both pages and either pag= may be used. 

The ASF Correction is addad algebraically to the 
time difference for the LORAN-C pair. Interpolation of this 
data will not necessarily improve the accuracy due to the 
method used to determine ASF Correctors [ DMAHTC, 198 1]. 
Since the correctors are computed in the azimuthal array and 
are based on the conductivity and distance over which the 
LORAN-C electromagnetic wave travels, =e wor. —_COnEeC TOR 
between the published ASP Correctors in the tables may not 
be the linear interpclated values. For example, eee Ae 
Mee=eccor for a distance of 500 m with an azimuth of 180° is 
equal to 1.5 psec. The ASF corrector for a distance of 500 
fee an azimuth or 1819 is equal t5 1.6 psec. The inter- 
polated value between 1.5 and 1.6 is 1.55. Dhice “a= ewA SF 
Corrector for the 180° 30' azimuth is 1.4 since the land 
distance for the same azimuth is less than the land distance 
for the 1809 and 181° azimuth. The LORAN-C signal passed 
Meet a harbor {Dansford, 1982]. 


ee TP TED DETERMINATION OF ASF CORRECTORS BY HYDRO fFfIELD 
EeR TIES 


One or the major problems encountered by paydrogrephic 
Seve y Operating units when using LORAN-C for posixzion 
Som=rO! 1s the determination of the AS? Correctors for the 
Baevey erea. The four surveys mentioned in the In*«roduction 
Menieaa= a=~empts to derermine the correctcrs by comparing 


moe LORAN-C rates te a second source. 





Bathymetric Surveys H-~9822 [NOAA 4H-9822, Hoole 1G. 
H-9823 [NOAA H-9823, 1979] Gulf of Alaska, compared the 
rates from an Internav LC-204 LORAN receiver to computed 
rates from a position obtained from a JMR-1 Satellites 
Navigation Receiver when available. Shore ties using radar 
ranges, visual bearings, and sextant angles in comparison ‘to 
LORAN-C rates were also made prior +9 and after each survey. 
The calibraticns of LORAN-C rates were based on the satel- 
ice positions only Since the positional computation of 
LORAN-C and JMR Doppler Satellite wer]? made on the WGS 1972 
datum whereas the land ties were basé¢d on tne NAD 1927 
datum. 

Bacehyme=ric survey SU-40-7-82 which extended along «he 
Washington, Oregon, and California Coasts used LORAN-C as 
navigational control. LORAN-C ‘time differences were 
compared with SATNAV positions. The raport did not indicate 

[NOA 


whether any correctors were applied A Shlp (Susvevor, 


mo 2). 

Finally, iy cesog =a pins Cc Survey S~K902-Wh-82, 
Reconnaissance Survey of Safety Fairways, Gulf of Mexico 
used LORAN-C as a positioning control systen. The posi- 
tioning unit was an LC-204 =receiver. LORAN-C rates were 
mmeut Via the HYDROPLOT Controller, a special purpose 
mpage ~-OutTput interrace which is the nucleus of the computer 
system nardware [Umbach, 1976}. Positions were computed and 
plotted by Program RK121, LORAN-C Real-Time HYDROPLOT 
meackus, 1980]. 

ASF Correctors for LORAN-C were achieved by visual cal:-~ 
bration uSing three point sextant fixes using charted oil 
mege as control in the vicinity of the survey area. A 
shree-point sextant fix is a convenient and accurate method 
for determining the position of a hydrographic survey 


vessel. sextants are used *oO measure two angles between 
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three objects of known geographic positon. The center 
object is common to both angles. The position of the obser- 
vers taking the angles is fixed by the intersection of three 
circular lines of position [Umbach, 1976]. 

These sextant angles were recorded and later tranferred 
onto their respective charts using a plastic three-arn 
Peo er ac.Or . A plastic three-arm protractor 15 transparent 
and made up of one fixed arm and ‘two movable arms which 
contain an etched line that is radial with the center of the 
Beotractor {[Umbach, 1976]. Sextant angles observed in pairs 
for a resection f1x with a common center mark mey be plotted 
directly by this instrument. When the three arms are placed 
at the angles observed and fitted sd as *o pass through the 
plotted positions of the observed stations on the field 
Sheet, the hole at the center of the three-azrm protractor is 
the fixed positon of the vessel [Ingham, 1975]. 

Partial correctors for each area surveyed were defined 
by comparing the observed rates and the determined crates 
me-.ed on the nautical chart. The partial correctors w 
applied via the HYDROPLOT Controller. However, even afte 
Mmeeeyeng «hese correctors, the plotted position still disa- 
greed with the ship's determined pestion with respect +o the 
Sel rigs. Ship's personnal attributed the discrepa 
one or more of the following: 

WekeCuracy Of the charted rigs, 

2) Weather effect on LORAN-C, 

3) Time of day, 

Deno pegaeaOr Dt Stgnal over land vata, 

5) Three-alm protractor accuracy, and 

SjmeeerOr in the conversion by the softwares 
of the LORAN-C rates to latitude and 
Pongsacuce (NOAA, 1982]. 
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There is an apparent need for a LORAN-C calibration 
routine aboard NOAA ships which provides the ASF Correéectors 
for program RK121, LORAN-c Real-Time HYDROPLOT. The routine 
should use the same geodetic distance computation found in 
RK121 and use the same datum as that of the nautical chart 
of the survey area. 

The above mentioned discrepancies illustrate the defi- 
Sm@encies in applying a single ASF Corrector to LORAN-C data. 
The accuracies for hydrography cannot be met using single 
cortectors because the errors are non-linear and systematic. 
They cannot be distributed like residuals ina traverse. 
Schnebele {1979} has already proven that single ASF area 
correctors +o LORAN-C positions do not meet the accuracy 
standards of the NOS Hydrographic Manual. 

Based on visual inspection 9f the DMAHTC LORAN-C 


Correction Tables, ASF Correctors should be updated every 
five minutes of Latitude or Longitude change. In Monterey 
bay, California, there 1s approximately 0.1 to 0.2 psec 
d:fference for every five minutes cf change, &@ potential 


meee Of 55 +0 110 meters. 
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III. EXPERIMENTAL PROCEDURE 
me FIELD PROCEDURES 
imeorder to Compare the use of differential LORAN-C with 
ASF multiple correctors, typical survey operations were 
planned hor tie southern portion of Montere Bay, 


California. This survey was conducted in conjunction with a 


evaluaticn of multiple lines of positon for 


1982}. 


on known geographic positions 


comparative 


selective positioning methods [ Anderson, FOwr —f=cEo— 


wave ranging systems were set 


Thesis vy ii 


Geographic Names and Positicns 


























i“. .  .... = oF ae tees aaa 

Microwave Systen ee coreh ass Osi econ | 

stations (NAD 1927) Datss used. | 

oe 4 (1964) 369 36' 23.44596" mame 3-5-1982) 3 

12704517 35.3328 1" 

| | 

ms Mon ECC. 869° 36° 04.7303)" June 3~5,1982 |] 

W21e S528) 35,930 50" i 

| 

Pe ceciver nee. 369 36' 32.49281" June 3—5, 19325) 

| Vane Soo e262 lal o 2 

| 

maussel Ecc. 36° 36' 18.25484" dime 3 — 5 ioe a 

1219 54" 11.49661" 

| 

Park (1931) 369 53 13.30600" June 6-7,1982 | 

1219 49' 46.74300" | 
Mulligan RM1 36° 44 56.49531" June 6-7,1982 

| VZOea7* 52,3 1090" | 

{Range 7 (1972) BO so Ol. 477387" June 6-7,1982 | 

| 1219 49t 98.58202" | 

Mussel (1932) Solas 1) 184 1o LOO" June 6-7,1982 j; 

1219 54* 1411.49661" 
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listed in Table VII. A series of tracklines were run in two 
separate areas as shown on Figure 3.1. To ensure that the 
microwave postioning systém was working properly, the equip- 
ment was calibrated over Known baselines of 1497.47 meters 
and 7877.31 meters at the beginning and end of the project. 
Trackline observations wer2 only mad= during the daytime in 
fair weather conditions s90 as to aliminate sky waves and 
weather changes that influence LORAN-C signal propagation 
@maracter:stics { Samaddar, Sis) ZANE The vessel used was the 
tee fcot R/V Acania which is operated by the Naval 
Postgraduate Schcol. 

The positioning equipment consisted of a Micrologic 
ML-1000 LORAN-C receiver (0.01 uséc resoiution) and 2 
Trisponder Microwave System orovided by Racal-DECCA Survey, 
Inc. The Trisponder Microwave System consisted of four DNT1 
Meme) 217C transponders, four DNT1 Model 21017 HP sector 
antennas with 879° by 5° beam widths, cne DNT1 DDMU (Digital 
Distance Measuring Unic), two Omni DVTI Model 21019 HP 
antennas, a Houston Instruments Model DP3-M2D/xRC3 plotter, 
and a Texas Instrumerts 743 terminal (Table VIII). The 
Mematecturer's published accuracy for the positioning equip- 
Meme 2S + 1 m for a single range [R8acal-DECCA Survey, 1981]. 
Anderson [1982] discusses ‘the accuracy of four lines of 
position. The four Decca Trisponder distances were recorded 
via a Texas Instruments 743 data terminal while the LORAN-C 
rates were manually logged. The data was acquired at one 
Minute intervals while tha ship maintained constant ccurse 
and speed. The recorded LORAN-C rates were 9940-W, X, ard ¥ 
Ot the West Coast chain. 

TO test the potential for calibrating the LORAN-C Systen 
uSing the Semi-Empirical Model, «the cerrection tables, and 
Multiple observed field correctors, the positions derived 

Ss 


= 
from the microwave system neaSurements were used to computs 
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Parone iV Le L 


Microwave Positioning Equipment 











pen oe Be ag a 
He Equipment szy | 
DNT1 Model 217C Transponders (Code 72R SIS | 
Code 74R 33720 | 
Code 76R 3321 
Code 78R Brae 2 | 
jDNT1 Model 2107 HP Sector Antenna 185 | 
| (279 by 5° Beam Widths) ee | 
| st | 
| 
DNT1 DDMV Model 540 426 | 
Omni, DNT1 Antenna Model 21019 HP 194 | 
{ 200 | 
Houston Instruments Plotter Model DP3-M2D/RC3 10722-10 | 
{Texas Instruments 743 Tarminal 34418 | 
Neen ne ee nn get cease ee eects ens cu “on “n-mes ems Sa Ga mn CSL lm Sa sD UDCA SSO 
expected LORAN-C time differences at each point. The 


difference or offset between these exvected time difference 


iF) 


and the cbserved values were computed for the tnaree methods. 
The mean offset, standard error, and drms values were also 


computed and compared. 


Pee MLCROWAVE SYSTEM POSITICNING 


The geographic position of the ship based on four lines 
of position was determined using a2 computer program callei 
MPsYLQ (Geographic Position by Least Squares) written by «he 
author (Se€e Appendix A). GPBYLQ contains subroutine LSQR 
(Least Squares), which is a least squares adjustment written 
Meeeeaut RK. Wolf, Ph.D. [1974] and revised by LCDR D. Leath 
imi |}. Geographic postions were converted *o X,Y (meners) 
Been in turn were converted t0 gsographic position via 
SubzoOutines GPTOXY and XYTOGP, respectively (Wallace, 1974}. 
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Subroutine GPTOXY and xXYTOGP are based on the Modified 
Transverse Mercator Grid which was centered in the survey 
area. 

The Modified Transverse Mercator (MTM) projection is 
used by the National Océan Survey and is Similar *o the 
projection used in the Universal Transverse Msrcator (UTM) 
system. The main differance is that in the MTM a Central 
Meridian is picked +hat is near the survey area instead of 
being fixed at a particular meridian [Wallace, 197 Wi. 
Central Meridian (CMER) , False Easting (ess), and 
Controlling Latitude (CLAT) are the three parametérs which 
define the MIM projection. CMER is the mean longitude 
computed using the maximum and minimum longitudes of the 
Buevey limits, FEST is the X-Coordinate that is assigned to 
the Central Meridian, and CLAT is the distance in neéters 
from the equator to seme reference latitude [Wailace, 1971]. 
The Céntral Meridian, False Easting, ande s€On21olling 
Latitude used for Monterey Bay, California referenced to NOS 
Maaet 18685 are: 

Gran = 121° 56" 00.0", 
woot = 20000. 0, 
CLAT 4050000.0, 


moO be consistent with the National Ocean Survey charts 


@e the areca, all computations were done relative to North 
American Datum (NAD) eZ SE GeOGEsOnsc posit ons. All 


programs were executed on an IBM 3033 computer located ar 
igi « Church Computer Center, Naval Postgraduate School, 


Monterey, California. 


eee LCORAN-C COMPUTATIONS 


The differences or offset between the observed and 
Computed LORAN-C rates using «he Sémi-Empirical TD Model, 
c 


a 
AoF LORAN-C Correcticn Tables, or the Multiple Observed 
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ieyebe 2X 


LORAN-C 9940 Chain Data 








Ae : Pema: 
; See Rac Position | 

| Staton (HaD® 1937) | 
|Master - Fallon, Nevada 399 33' 07.03"N { 
1189 49* 52.23"W | 

Slave - eer oe Washington 479 Q3' 48.82"N | 
l 9946-w 1190 44t 34:78"W | 
| slave - Middletown, 389 46* 57.49"N | 
! Caletioumaa 99 40=x% 1229 29* 4O.O4"W | 
|Slave - See Ce” Nevada 359 19* 18.32"N | 
| 9940- 114° 48* 13.95'WH { 
eer 


Correctors were compared to the offsets between the coserved 
TD rates and the computed rates fer which only the s2awater 
Secondary Factors (SF) were applied. The comparison of the 
offsets between the four mathods illustrates the improvement 


maePOsSLtional accuracy after applying ASF Correctots. 
1. Seawater secondary Factors (SF) 


Time differences using only the s¢awater Secondary 
Meetors for each of the geographic positicns were computed 
Hoang program LORAN written by the author (Appendix C). 
seawater Secondary Factors (SF) wer2 computed uSing formula 
fro and the coefficents found in Table VI. Ad TO4 
distances in meters were determined using subroutine INVER1. 
Mier! 25 a geodetic inverse routine using T. Vincenty's 
modified Rainsford's method with Helmert's elliptical <erms, 
programed by LCDR L. efeifer, NOAA [1975]. Subrecutine 
INVER1 is accurate to 0.0901 m halfway around the world 
Merciter, 1982]. All distances were converted to microse- 
Conds using 299.792458 m/sec. Time differences (TD) w 


—- 


— he 
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computed from equation 2.5. North American Datun 1927 
geographic positions were used for all computations. See 
Table IX for the positions of LORAN-C 9940 transmitters 
mazvordan, 1979]. 

To ensure that subroutine INVER1 was functicning 
properly the distances between the master and secondary 
stations were compared to the NOS published baseline 
distance [{Riordan, 1979]. The published distances and the 
results from routine INVER1 are listed below. 


Published Baseline Computed Bassline 
Distance (nm) Deseat ea ga ewe RA | (a) 
9940-W 837,777. 0929 837,777.115 
9940-X 327,886. 3720 327,886. 216 
9940-Y 589,298. 5712 589,298.589 


The difference between the published and computed basélines 
meang]ed from 0.02 to 0.06 n. 

Differences or offsets (x. ) were obtained by 
Bemoetracting the observed LORAN-C rates from the computed 
values from the various nethods. The mean difference or 
@eeset (x) and standard deviation (s) in microseconds for 
each rate were determined using equations 3.1 and 3.2 
Mmeoennacott, 1935]: 


Mean offset (psec): Rea eK. (35 1) 


=a 


Wie mewex es = OL gn al obse=vaclon in usec, 
1 
nm = number of observations; 
tandard 1 o 
Soe aloha Say > exes © (2 - X)¢. (3.2) 
(usec) Rel. 2 
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The mean offset in microseconds can be converted to meters 


using equation 3.3 [Heinzen, 1977]: 


Mean offset xX wW 
(meters): = 2s a (333) 
Sin a 
where w = distance corresponding to oné microsecond 
on *h2 baseline 
= 1¢97a90229 Meters [| Bigelow, 1963]-; 
= ene-nal = sens angle between the | 
radius vectors from the position to 


mamation 3. 


{ Heinze 


drms= 


Mae, correlation coefficient (9p) is a resu 


Beaticn 


prea ton 


the master and secondary stations. 


4 was used t> compute drms values in meters 
Mee 19172 





Ss w)e S w)* 29cos(a +a Swsw 
1 ( 1 2 ( 1 2? 1 Z 
oo ----- -- ---- torn ere ter meer ere rrr (3. 4) 
Sin( + ) 2 2 Sine y (Sana } 
1 2. Sina sin a 1 
1 2 
where: Oe CoOspelatmeon cosfftberent = 0.33, 


S and a ar2 aS above with the subscript 
denoting the appropriate 
secondary station. 


f 


}- 
cf 


of the secondar 


ition with the master 


qt 


having a common line of »po 


" 


‘ Although often ignored, various authors assign 


values ranging from 0.33 to 0.40. Bigelow {1963] chocses 


9 = Q. 
- 
reduce 


Se 
semi-Empirical ID Grid 
Powedcccem@ine 2£ the Semi-Enbirical TD Grid would 
€ither the drms value or the offset between the 
observed and expected TD rates, S-egran LOPLC (Lins. <of 


Positio 
Progran 
Water, 


=he meee 


mn - LORAN-C) was written by the author (Appendix B). 
LOPLC computes the distance over land, distance over 
PiomeotralGsSscance, and ths azimuth from north for 


ansmitting station using Subroutine INVER1. The 


4 
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land-sea distances were computed by selecting 4 point which 
was located along the coast of Monterey Bay. Subroutine 
SHORPT (Shore Point) interpolated a geodetic position from 
Memageogetic points which outlined Monterey Bay, California. 
The land/sea Secondary Factor was d2termined in Subroutine 
SECFAC which is found in program LOPLC uSing 2quations 2.8 
eerough 2.12. All computations were based on NAD 1927 
geodetic datum. 

Tests were made on pregram LOPLC using data found in 
mae Semi-Empirical TD Grid article [The Analytic Science 
Corporation, 1979] using the WGS 1972 datum. Stat LOnmnASC 
55 located at latitude 349 34 18.3" N and longitude 120° 
39° O25 WW, was selscted from The Analytic Science 
@erporation article. It was one of the stations used to 
calibrate the coefficients for the Sami-Empirical TD equa- 
tions discussed earlier. The only distances listed for TASC 
55 were the individual distances over land and over water. 
The total distance between the transmitters and TASC 55 was 
computed by adding the land and sea distances. Tha 
BeelOWLNgG Station-tc-site path segment lengths are listed 
meme TASC 55; 


TASS 35 INVER 
Land 32a Tozal Total 
Distance Distance Distance Diszance 
Station (kK) (km) (kK a) (km) 
Master 540.730 Soo ue 5 765 270 576.083 
X-Secondary 369.248 126.096 495.344 495.178 
Y-Secendary 525.659 15.832 541.491 Sues OS 


Using the land distances, sea distancé, and soztal distances 
femtinoC 55, and the computed azimth from subroutine INVER1, 
program LOPLC produced offsets of ).86 prsec and 1.13 psec 
ror 9940-X and 9940-¥ respectively when compared +o the 


expected time differences at TASC 55. 
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An attempt was made to determine 1f this discrepency 
was due t¢0O program LOPLC. A comparison of the total 
distances in the report between the transmitters and TASC 55 
were compared to the total distances computed by INVER1. 
The difference in distance between Subroutine INVER1 and thea 
total distances from the Semi-Empirical TD Grid acticle 
Genged from 160 - 200 m corresponding to 0.5 to 0.65 usec. 
Maes results in time difference errors of 0.04 usec for the 
mee0-X rate and a 0.10 psec for the 9940-Y rate, implying 
faae the offsets, 0.86 psec and 1.13 psec, are caused in 
part by the method in which the azimuth from north is d¢éter- 
mined. Unfortunately, azimuth data from TASC 55 was not 
presented in the article. The Analytic Science Corporation 
has been contacted on numerous occasions in an attempt to 
ascertain their method of determining distance and azimuth. 
As of this date there has been no response. 

Nevertheless, the data from Monterey Bay was 
[iemltzed in Program LOPLC to obtain results that could be 
compared to that obtained by the] other methods. BE -chis 
method is accurate enough, ASF Correctors could be deter- 
Mined via computer for each individual position without 
uSing tables or field determined correctcrs. The mean 
Meeset (xX) and standard error (s) in microseconds between 
the observed and calculated rates were computed using equa- 
mons 3.1 and 3.2, respectively. The mean offset (1) and 
drms in meters were computed using equation 3.3 and 3.4, 


respectively. 


Pee calculated Table ASF Correctors 


=a. —_ —_— ——_ =) oe eee ee ee ee eee ee ee 


The offset between the cbserved LORAN-C races and 
the expected values with applied ASF Correctors from the 
EORAN-C Correction Tables and «he seawater Secondary Factor 


(equation 2.13) were determined using program LORTAB which 





was written by the author (Appendix D). Memon COrLeCtor 
for each position was salected using subroutine TABLE. 
Subroutine TABLE, which is found in program LORTAB deter- 
mines an ASF Corrector for each data point based on its 
geodetic position. The ASF Correctors used in subroutine 
TABLE (see Figures 2.5, 2.6, and 2.7) are located between 
latitude 369 35' N and 36° 55'N and longitude 1219 50' W and 
m22° OO dW. The difference in sign between the ASF 
Merrectc=s in subreutine TABLE and those found in <he 
LORAN-C Correction Table is due to the difference in their 
pep lication. ASF Correctors from the tables are applied to 
observed rates while ASF Correctors from Subroutine TABLE 
are applied to the calculated LORAN-C rates. Negative ASF 
Merrectors from the LORAN-C Correction Table were applied to 
the calculated «ime differences to be consistsn= with the 
application of Secondary Factors +> the computed oprimary 
phase delay. 

As before equations 3.1 through 3.4 were used ¢t 


O 


th 


@empute che the mean offset, standard error, and drms. iz 
Mare application of LORAN-C Correction Tables is accurate 
enough, i* precludes the need to determine ASF Correctors in 
the field. 


4. Observed ASF Correctors 


Observed ASF Correctors were determined using 
Program ASFSEL (ASF Selection) which was written by the 
Suenor (Appendix E), Schnebele's prior data [1979], and «he 
June 1982 data. Program ASFSEL (ASF Sélsction) was written 
by the author. Dice lEojsaenmeconpuces che ASF Correc*tors by 
Subtracting the observed LORAN-C rates from the expected 
values. Only the seawater Secondary Factors from equation 


2-13 have b2en applied to the calculated time ditferences,. 





The mean ASF Correctors for the LORAN-C rates were deter- 
mined at every minute of latitude and longitude between 369 
50' N and 369 35' N and 1229 QO4'W and 1229 49' We. See 
Appendix E for mean ASF Correctors at one minute intervals. 
The ASF Correctors were then selected and assigned to 


subroutine TABLE in Program LORTAB at five minute latitudes 


and longitude intervals. All values were entered to the 
nearest hundredth of a microsecond. Equations 3.1 through 
3.4 were used as before for computations. This determina- 


meen Gr ASF Correctors in the field, if accurat2 enough, may 
allow the use of LORAN-C as a positioning system for hydro- 


graphic surveying in the future. 


oe 








IV. RESULTS 

A total of 620 time differences (TD) and geographic 
positions based on four lines of position were recorded in 
the southern portion of Montery Bay in order to compare the 
use of differential LORAN~C with ASF Multiple Correctors. 
The data was divided into four sets. The first data set is 
Bennebele's [1979] prior data consisting of a total of 130 
data points collected on two separate days, Jun2 12 and July 
=, 1979 between latitude 369 38' N and 369 47" N= and 
longitude 1219 49* W and 71229 OQ2' W. The recorded LORAN-C 
time difference rates were 9940-Y and 9940-W. 

mMhe seccnd data set consists 9f 193 time dirferences 
collected on June 7, 1982. The recorded LORAN-C rates were 
also 9940-Y and 9940-4. Data set Number 2 is located 
between latitude 369 4O' N and 369 45' N and longitude 1219 
54* W and 1229 OO'WN. Sssperratese 441) Fos the .locar2on of 
test areas in Monterey Bay, California. 

Data set Number 3 is located in the same area as data 
set Number 2, between latitude 36° Uo? N and 369 45' W and 
longitude 1219 W 54" and 1229 OO W. Mis et, CONS. Scing 
Beeeizo points with recorded rates 9940-K and 9940-Y, was 
recorded on June 6, 1982. 

The final set, data set Number 4& was recorded between 
Jun¢ 3 and ,June 5, 1982. femGouet2ns. 169 polnts located 
between latitude 369 36 N and 369 39" N and longitude 121° 
Boe W and 1219 58* W. Pie recon tegmra.es were 9940 —xX and 
9940-Y which are the same as data set *hree. The data for 
rates 9940-X and 9940-Y was kept separate so 3s to determines 
if there was a significant difference between «he ocffshors 
(data sst three) and inshore (data se¢* four) drms values due 


tO Dhase recovery (see Chapter Two). 
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TABLE X 


Data Set Parameter 

















[ Set Collection TD Number of Area Limits | 
) 224258 Date Rates Data Points Lat Lon | 
{ 1 June 12, 1979 9 940-Y 130 36/38 1594031 
| July Zone 1979 9 94U0-WN 36/47 aimed! 
n 2 June 7, 1982 9940-Y 193 36/40 121/54] 
9940-W 36/45 122/00} 
3 June 6, 1982 9940-X 128 36/40 121/54 | 
9940-Y 36/45 pec 

4 June 5, 1982 9940~ xX 169 36/36 1217 53) 
9940-Y 36/39 121/58| 

eee eee eee 


Table X provides a convenient breakdown of the parame- 
ters for each of the data sets. The table consists of the 
data set number, the date the data was collected, +he 
LORAN-C time difference rates, the number of data points, 
and the area limits in latitude and longitude. It detines 
the parameters for the four data sets of Tables XI (Seawater 
Seeeomdary Factor Errors), iE e(S]2mi-Empirical TD Grid 
Semprec-ion Errors), XIII (Multiple LORAN-C Correction Table 
meeoms), and XIV (Multiple Observed Corztection Errors). All 
basic data and computations are kept on file with the NOAA 
hydrography instructor at the Naval Postgraduate School, 
Memcerey, California. (Individuals seeking this information 


Should contact th2 Oceanography Department.) 


A. SEAWATER SECONDARY FACTCRS 


LORAN-C rates were computed using only the ssawater 
Beeendaty Factor (equation 2.13). The LORAN-C computed time 
differences are basically uncorrected rates Since no 


Additional Secondary Factors were applied. Offsets between 
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TABLE XI 


Seawater Secondary Factor Errors 


————-—SO eine ae one on aeeee = oa 
| Mean Mean Standard 
| Offset Offset faa Of drms { 
Rates (x-s ec) (1-m) (S-ps¢c) (m) 
| 
| a) Set 1 (Schnebele's data, 130 data points) { 
Y -0.508 -~148.6 36 ).088 110.4 
| W -1.241 -6 83.693 0.134 | 
by set 2 (June 7, 1982, 193 data points) \ 
Y -0.526 Sue 9S Ome 7 Ones 
W -1.283 -701.947 0.124 | 
| eset 3 (June 6, 1982, 128 data points) | 
X 1.565 5 35.449 0.059 38.6 
| Y -0.550 -161.404 0.083 | 
d) set & (June 3-5, 1982, 169 data Bae | 
X 1.582 §50.089 0.05 68.1 
Y -0.817 ~-237.109 O23 7 
= peau hn | 


the observed and calculated rates using only “the seawater 
Secondary Factor were computed to illustrate =he in vement 


pro 
mepOsS2-i0on aftexz applying the Semi-Empirical TD Model, ASF 


Deaen-C Correction Tables, or the Multiple Observed ASP 
Correctors. The mean offset in microseconds and neters, 
standard error, and the drms are found in Table iXI. 


Examples of offsets for several data points are listed in 
Appendix 8 after program LORAN. 

It was stated earlier in Chapter One that Schnebels 
obtained a2 66 m 1 drms using Y and W crates. The drms of 606 
Mwas obtained using 48 data points which were located 10 kn 
or more offshore. The 110.4 n 1 drus for the Y and W rates 
in Table XI is a result 9f combining the 130 inshore and 
offshore positions. The increase from 66.0 m 1 drms for the 
MesenOore positions to 110.4 m for combined offshore and 


inshore positions indicates that the application of a single 
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ASF Corrector is dependant on the size of the area. The 48 
data pceints were located between latitude 36° 41" N and 369 
mer N and longitude 1279 55 Wand 1229 O2* ¥W. 

Schnebele also obtained the 66 m 1 drms by adjusting the 
LORAN-C observed time differences which were skewed due to 
the ship's motion and the five second averaging interval of 
the LORAN-C receiver. These caused the observed TDs +o be 
several seconds old in comparison to the microwave sys*en 
measurements [Schnebele, 1979]. Due to the large amounts of 
data from the June 1982 survey operations, no deskéewing was 
done. 

The offset, standard error, and drms for data se+ Number 
2 (June 7, 1982) compares well to data set Number 1 
(Schnebele's data - all). ISO pweaeS Soeoan | drms for the 
X-Y rates is between Schnebele's 42.0 m 1 drms prediction 
{[Schnebele, 1979] and Nelson's findings of 30 m 1 drms in 
Mmaerrancisco Bay [General Electric I0o., 1979}. The largs 
mens Of 68.1 m for data set Number 4 (X-Y¥ rates) is probably 
due to the phase recovery of the electromagnetic wave from 
the 9940-¥ transmitter located in Searchlight, Nevada. 
femal iy, if ASF Correctors are not applied, drms values 
Memoged trom 38.6 m for 9940 X-¥Y +0 101.4% m for 9940 Y-W for 
Gata sets 2 and 3. The large offsets for all data sets 
indicate a systematic error, ranging from 150 m to 700 qn, 
Which precludes the use 9f ‘this method for hydrographic 


surveying. 


Pee SemMi-EMPIRICAL TD GRID 


To determine if the Semi-Empirical ID Grid would reduce the 
@eeset and drms for hydrographic surveying, program LOPLC 
S$ 1 through 4. The 


Semi-Empirical TD Grid apoliss a Secondary Factor and an 


(Appendix C) Was applied to data set 


Additional Secondary Factor +o the primary phase delay based 


oT 





on the distance overland, «he distancé over water, and the 
total distance using mean land and water conductivities. 
The mean offset in meters and microseconds, standard error, 
and drms are listed in Table XII. Examples of offsets for 


several data points can be found in Appendix C after proqram 


TABLE XII 


Semi~Empirical TD Grid Correction Errors 











= nn aren ar =) 6. ear eve agen ep toe ee 7 
| Mean Mcan Standard | 
| Offset Offset Ear dx-ms | 
Rates (xX-psec) (l- m) (S-~nséc) (m) | 

| 

a) Se+ 1 (Schnebele's data, 130 data points) | 
Y eos 330.672 0.086 ssa ae 
| 1.745 650.278 Onto So | 
may set 2 (June 7, 1982 data, 193 data oe 
Y eos 324.975 noua 1702 7 atl 
| W eS 603.926 Om tZo | 
ieoy Set 3 (June 6, 1982 data, 128 data ee { 
ie on 386. 940 .061 3980 a 
| Y 0.589 W239 23 0.084 
| ad) Sex 4 (June 3-5, 1982 data, 169 data Reon >? | 
X 1.050 365. 113 0.052 68.2 | 
i ie 0.341 Soe ca 0.189 
—- eens cies aerailecmeeiees teins cic ets daiainis eves Sia gece enemas a aoc aaa 
HOP LC. Earlier, Program LOPLC had been tested with data 


Memnad 25 the article by The Analytic Science Corporatio 


3 


mo 79 j. Program LOPLC produced offsets of 0.86 msec and 
1.33 psec for the time differences from rates 9940-X and 
moe0-Y, respectively. AS previously stated, «he large 
offsets may be due to the method by which the distance and 
azimuth were computed. Th2 drms, for the data from Monterey 
Bay, Obtained with the S2mi-Empirical TD Grid was similar to 
the drms errors for séeawatar Secondary Factcr Model. is 

re 


Maght be an indication that Program LOPLC is correcz but the 
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bias needs to be adjusted to reduce the large offset between 
the observed and calculated [TD rates. The mean offset 
ranged from 99 to 630 m. Bias reduction could be achieved by 
applying land-sea data for the entire West Coast +o the 
Semi-~Empirical model. Again, the existence of large offset 
values precludes the use of this method fer hydrographic 


surveying. 


ee LABLE ASF CORRECTORS 


To determine if the multiple ASF corrections from the 
LORAN-C Correction Table would reduce the offset and drms to 
meet the NOS accuracy standards, program LORTAB (Appendix D) 


was applied to data sets 1 through 4. RSP COrece Ors tron 


TELE XEae 


: Miele obeeLORAN-C Corzeetion Table Errors 








a. a ise a ee 

| Mean Mean Standard | 

| Offset Offset nike or drmms i 

! Rates (x=psec) (i-m) (S-nsec) (m) | 

| 

Meaye set 1 (Schnebele's data, 130 data points) 

| ¥, =0.229 -67. 068 0.104 123.6 | 

W 0.169 92.540 0.160 | 

| Dye set 2 (June 7, 1982 data, 192 data points) 

| { =0.309 -96.614 0.085 116.9 | 

W Os hes 67.189 0.144 | 

| Syescet 3 (June 6, 1982 data, 128 data Soe), | 
X 0.559 194.134 .073 47.3 

: a -9.319 Sosa o On 107 | 

ayeeces 4 (June 3-5, 1982 data, 169 data points) | 

NemeotrCOnDpcetz ons Listed in Table for the south- | 

east end of Montersy Bay next to the shore iine. | 











pe Hirai elegy dolled I 





P= —eEmEomAN=-C Correcticn Tables are determined fron fiel 


fas yh 


Seer va-20n Of land conductivities by «he U.S. Coas= Guar 
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Calibration Team [Marine Science Department, 1982 ]. The 
results are listed in Table XIII. 

When compared to the Seawater Secondary Factor Error in 
Table XI, the drms values using LORAN-C Correction Tables 
were increased slightly while the offsets were reduced 
substanziallv. The offsets ranged from 65 - 200 nm. Since 
mae drms ranged from 47.3 =o 116.9 m for rates 9940 X-Y and 
9940 Y-W respectively, the application of ASF Corrector from 
the tables does not meet the NOS accuracy standard. Again, 
large offsets and the increase in drms precludes the use of 


mmas method for hydrographic surveying. 


tee MULTIPLE OBSERVED CORRECTORS 


To determine if multiple observed correctors would 
diminish the offse+ and drms values, mean ASF Correctors 
were selected at one minute latitude and longitude intervals 
uSing Program ASFSEL. One minute ASF Correctors are shown 
at the end of Program ASFSEL in Appendix E. From the one 
minute grid, mean ASF Correctors were seiected and entered 
into Subroutine TABLE at flve minut2 Latitude and longitude 
intervals in Program LORTAS3. The following is an example of 
Maeeg940-X ASF Correctors at five minute latitude and longi- 


tude intervals for the program: 


122/00 /00 .0 UZ 7200-2 0 121/50/00.0 
36/50/00.0 at 2 = eos 
36/45 700.0 ae 2 Sle Sole 
36/40/00 .0 sso = eo) 


This is the same format used in the LORAN-C Corzection 
Tables. Program LORTAB was applied to the four data sets. 
The error results are listed in Tables XIV. 

The drms value obtained with multiple observed correc- 


sors were all reduced when compared to the adrms for the 
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Woon kV 


Multiple Observed Correction Errors 








i - ee 
Mean Mean Standard | 
| Offset Offset pel Or drms | 
| Rates (X-psec) i=) (s-psec) (m) . 
| a) Set 1 (Schnebele's data, 130 data poaeae 
y 0208 3 9-305 -076 B75 | 
i W 0.052 28.604 C23 | 
meeset 2 (june 7/7, 1982, 193 data points) { 
Y 0.028 8.318 0.073 89.3 | 
W ~0.016 ~8.944 0.116 { 
eyescet 3 (June 6, 1982, 128 data poin<s) 
X -0.034 ~ 11.490 0.055 38.3 | 
i y 0.010 22913 C2096 | 
meoet & {June 3.5, 1982, 169 data eae { 
{ X 02007 0.383 0.05 67.5 
Zz v ~9.008 =25225 Oueiey 





seawater Secondary Corrector. The most impressive reduction 
in drms was within Schnebele's data which covered an area of 
seven minutes of latitude and 12 minutes of longitude. The 
drms for seawater Secondary Correctors was 110.4 m whereas 


the drms for the same data using multiple observed corrzec- 


tors was 87.6 nO. This is a smaller drms than that of the 
June 7, 1982 data (data set Number 2) which was obtained 
three years later. ieamapoears that 37.6 m 1 drns is nearly 


Or Vene so iO" y= h 


the minimum error that can be obtained 
rates in Monterey Bay artes applying multiple observed ASF 


- 66 


Beefectcrs at tive minute latitude and longitude intervals. 


For the LORAN-C rates 9940 X-Y, a 38.3 mn 1 dzrms was obtained 


ho 


for same five minute area covered by data set Number (June 


foe 1982). LORAN-C rates 9940 X-Y were not obtained for tha 


ct 


Same size area covered by data set Number 1 (seven minute 


i) 


Of latitude and <welve ninutes cof longitude - 9940 Y-W) dua 
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mo the restriction of ship time and the length of time the 
Racal-DECCA Trisponder electronic equipment had been loaned. 

The drms for the June 3-5, 1982 inshore data was only 
reduced to 67.5 m from 68.1 m for seawater Secondary 
iGercectors. The small change in error at the coast is 
probably a result of the erratic behavior of phase recovery 
discussed earlier in Chapter Two. The drms value couid 
presumably be reduced if the correctors were applied at one 
Minute intervals. This would be 32 very costly method of 


calibra*ing Loran-C for hydrographic surveying. 
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V. CONCLUSIONS 

Tt was noted ir Chapter One that the smallest scale 
routinely used for coastal surveys is 1;80,000. This yields 
an allowable error of 49 mn 1 drms with no systematic errors. 
This paper determined whether or not applying multiple 
Meet ional Secondary Factors (ASF) CORE eCCOES ro S LORAN=C 
lines of position would reduce the drms sufficiently to meen 
the accuracy standards set by the National Ocean Survey. 

Three methods of applying multiple ASF Correctors were 
tested. The first approach computes the *ime differenca 
based on a Semi-Empirical [ID Grid. The Semi-Empirical Model 
produced large offsets in the 9940-W, 9940-X, and 9940-Y 
ime differences. The offsets ranged from 99 to 630 a. The 
drms for 9940 X-Y combination was 39.4 m and the drms for 
9940 W-Y combination was 102.7 2. 


The second method applies ASrf Correctors found in the 
BeeHTC LORAN-C Correction Tables =o. LORAN]C. tines? Of 
position. The application of the tables reduced the offset 
in the LORAN-C time differences. The offsets wers between 


67 and 191 m. The drms was increased to 47.3 n for the 9940 
mar CcCOmbi=nation and 116.9 m for the 9940 W-Y¥ pair. 


Mhe final and most accurat2 me2thod applies multiple 
observed ASF GCOEESEr9ES oat Five: Minute latitude and 
Memagitude intervals to LORAN-C lines of position. This 
me=nod again reduced the offset in the time difference. 
This offsecr was between 3 and 12 m “for the 9940 X-¥ 
eemoination. Part of «the offsat may have been a result of 
the microwave positioning systen. Reference is made <*0 

p : 


Anderson's [1982] paper (inpre 


Meco ple lines of position. 
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The drms values were also reduced to 38.3 m for the 9940 
K-Y rates and 89.3 m for the 9940 W-¥ combination. The 38.3 
m1 drms can be decreased by improving the sampling time for 
LORAN-C receivers. Nelson obtained 30 m 1 drms for the 9940 
X-Y cates with special LORAN-C equipment used in San 
Francisco Bay, Gaiitornia {General slectric Co., 1979 ]. 
Improving the sampling time for LORAN-C receivers used as 
positioning equipment for hydrographic surveys should be 
investigated. 

heel Ggnams Values of 38.3 m with the possibility of 
obtaining 30.0 m1 drms and offsets ranging from 3 to 12 on, 
it may be possible to use LORAN-C for hydrographic surveys 
at scales of 1:80,000 or less using multiple observed ASF 
Correctcrs. The use of DMAHTC LORAN-C Correction Tables 
should not be ignored. After updating these ASF Correctors 
with observed data, the LORAN-C Correction Tables nay allow 
LORAN-C =o be used as a positioning system for hydrographic 


surveys. 
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PROGRAM GPBYLQ 


GENERAL PROGRAM FOR DETERMINING GP FROM KNOWN STATION 
BOoeETONS ANDTHE PISTANCES FROM, THEM OSING LEAST 
SOUNKEo.s Ge TO XY°AND XY £0 GP ARE DETERMINED BY 
SUenOuleN bs GeroxY AND XYTOGP /WHICH ARE BASED ON THE 
DO Dir eoUmltANoveRos MSRCATOR PROJECTION (MTN). 


PROGRAMMED BY GERALD E. WHEATON, LT. NOAA 


Meno tl SOUARES A DYUSTMENT BY PAUL R. WOLF, PH.D. AND 
BeemoeD BY Ds DEATH, LCDR 


PROGRAM INPUT VARIABLE NAMES 

Mii ANY JOB IDENTIFICATION NAMES OR NUMBERS 

M AND N = THE NUMBER OF EQUATIONS (M) AND UNKNOWNS (N) 
ROMANO YO >= BEST QUES TIMATE OF THE POSTION 

Sai, )) = x COORDINATE OF KNOW STATION 

SlAt(l - 2) ¥Y COORDINATE OF KNOW STATION 

anys) = DESTANCE FROM KNOW STATION 

Meee) = LHE COEFFICIENT MATRIX 

Eee pe) =] THE CONSTANT MATRIX 

OLL(I,J) = THE WEIGHT MATRIX (WEIGHTS ARE ENTERED AS 
ooo eeeSOLUTLON TS EQUALLY WEILGHTED) 


IMPLICIT REAL*8 (A-H, 0-2) 

DIMENSION XCORD (30), YCORD(30) ,¥1T (30) , ISNO(30) 

CCMMON /ISTAT/ STA (30,3) ,PXY(30) 

COMMON /LSOX/ A (30,30) ,EL (30,1) ,QLL (30,30) ,AT(30,30), 
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ome OV On Kt 30,350) , AQL(30,10), X(30,1), V¥{30,1), 
2VAR(30), TITLE(80) 


c 
S READ AND WRITE OUTPUT TITLE 
c 
WRITE (6,509) 
READ (5,710) TITLE 
WRITE (6,710) TITLE 
c 
€ DEFINE NUMBER OF SISNALS (NOT GREATER THAN 30) AND 
c NUMBER OF DATA SETS. 
C 
READ (5,502) NSIG,NDATA 
502 FORMAT(I5,110) 
C 
C DEFINE THE CENTRAL MERIDIAN (CMER), FALSE EASTING 
c (PEST), AND CENTRAL LATITUDE (CLAT). 
c CMER IS EXPRESSED IN DEGREES, MINUTES, AND SECONDS. 
c FEST IS THE X-COORDINATE THAT IS ASSIGNED TO THE 
Cc CENTRAL MERIDIAN AND IS EXPRESSED IN METERS. 
c CLAT IS DEFINED AS THE CONTROLLING LATITUDE. 
c IT IS USED TO REFERENCE THE Y-COORDINATES AND 
C IS EXPRESSED IN METERS. 
c 


READ (5,503) ILONC,ILMINC, RLSECC,FEST,CLAT 
es) FORMAT(1X,13,1X%,12,1X ,F8.5,1X,£7.1,F10.1) 
CMER = ((IABS(ILONC) * 60 + ILMINC) * 60) + RLSECC 


DEFINE VARIABLE FOR: 
NUMBER OF EQUATIONS (M) 
NUMBER OF UNKNOWNS (N) 
IPAGE = NUMZER OF LINES PER PAGE. 


qaaanaana 
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NYA A A w 


800 


805 


NULL WEIGHTS 


DO 3. IM=1,30,1 
DO 3. JM=1,30,1 
OLL(IM,IJM) = 0.0 


Ponape olm.lON NOMBERS ,POSITION, AND WEIGHTS. 
menvenie ©rOSuTLOMS (GP) TO SECONDS AND THEN TO XY. 


DO 12 J=1,NSIG,1 
Be (3, 600) LoSNO(d) ,TLAT, TL MIN, RSEC, JLON,JMIN,SSEC,NWT (J) 
PieemetkK, Lope ls,r9 .5,15,13,F9.5,F5.1) 


RMAST = (({IABS(ILAT) * 60 + IMIN) * 60) + RSZ=C 
RMASTL = ((IABS(JLON) * 60 + JMIN) * 60) + SSEC 


CALL GPTOXY (RMAST, RMASTL, XMETER, YMETER, FEST, CLAT, CMER) 
XCORD (3) XMETER 
YCORD (J) YMETER 


Pew ob e(STATION NUMBERS AND THE DISTANCES) 


DO 40 JCOUNT=1,NDATA,1 

peo, COS) LF,IiS,IT, 14,NREC 

PO RMAT (515) 

Bometo I{COUNT=1,NSIG, 1 

meat oNO(LCCUNT) ~EQ. IF) IF=rCOUNT 
PeeeononeeCOUNT) 220. IS) 2fS=ICOUNT 
Par oNO{LCOUNT) «80. IT) IT=LICOUNT 
ir (ISNO(ICOUNT) .2Q. I4) I4=IcouNT 
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507 


20 


SA, 1) 


MEOR Dir) 


STA(1,2) YCORD (IF) 
Omi 1} WT (IF) 
STA(2,1) XCORD(IS) 
STA(2,2) YCORD(IS) 
Orne 22) aT (IS) 
STA(3,1) X¥CORD(IT) 
STA(3,2) YCORD(IT) 
QLL(3,3) WT (IT) 
STA(4,1) XCORD(I4) 
STA(4,2) YCORD(TI4) 
QLL(4,4) WT (r4) 


READ THE DISTANCE RECORD AND LORAN RATE. 


DO 38 KCOUNT=1, NREC,1 
READ (5,507) STA (1,3) ,STA(2,3) ,STA (3,3) ,STA (4,3), 


1 BATE1, RATE2 


PORMAT(6F10.1) 


Doren eNe BRST GUESS COORDINATES XO AND YO 
Weta suUBROUTZNE GUESS. 


CALL GUESS (X0O,YO) 
IGUMP = 0 
COMPUTE FXY, 


A AND L MATRIX 


BOs25 [-1,% 
Pet) == DSORT (DABS(({ AO-STA (LI, 1)) **¥2 
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a) 


Ya aA | 


30 


Bo 


1 (YO-STA(I, 2) )**2) ) 

A(I,1) = (XO-STA(I,1)) / (PXY(I)) 
A(I,2) = (YO-STA(I,2)) / (FXY(I)) 
EL(I,1) = STA(I,3) ~- FXY(I) 

CALL SUBROUTINE LSOR 

CALL LSOR (",¥) 


COMPUTE THE NEW QUESSTIMATE FOR XO AND YO 


xO = XO + X(1,1) 
YO = YO + X(2,1) 


EXIT IF STANDARD ARE MET USING IJUMP OR 
AOAND YO CUT OFF 


Momo (tk tt, 1p) es Lhe le OO 2AND. DABS (X(2,1)) .LE. 1.090) 


fae GO TO 35 


fmeqtgumPe <EQ. 10) GO TO 35 
GO TO 20 


COMPUTE ERROR ELLIPSE 


CONTINUE 
CALL ELIPSE(SU, SV) 


SONVERT Xi TO GP 

melee XY LOGE (xO, 10, SEC LAT, SECLON, FEST, CLAT,CMER) 
Mlb LOpms (SECLAT, EDEGP,IMINP,RSECP) 

GCaALL TODNS (SECLON, JDEGP,JMINP,SSECP}) 


Poeun DeCONMENUE WITHA NEXT SET CF OBSERVATIONS 
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Web EMOgovul) IDEGE  PNeINEORSECP, JDEGP, JMUINP,SSECP, 
7 RALE |,Z,RALEZ 


C WREDE (6,505) IDEGPSIUINP, RSECP, JDEGP,JMINP,SSEC?P, 
asUu, SV 

505 rama (L483, F6.2,15, 13,F6.2,2F9. 2) 

ie Pe(PPAGHes tO. 50) WRITE (6 ,509) 
IF(IPAGE .£Q. 50) IPAGE = QO 

38 IPAGE = IPAGE + 1 

S 

C FORMAT STATEMENTS 

c 


710  FORMAT(80A1) 
509 FORMAT(1H1) 


40 CONTINUE 
STOP 
END 
Css 2555552 S2525 555525 252555 2SS5 SSeS S55 S25 555555555 55255555525 


nee eo ee ee eee 
IMPLICIT REAL¥*8 (A-H,0-Z) 
COMMON /LSQX/ A (30,30) ,£L (30,1) ,QLL (30,30) , AT(30, 30), 
mies, 30), OF n0s0ro0) , 200(30,10), X(30,1), V¥(30,1), 
2VAR(30), TITLE(80) 
COMPUTE A TRANSPOSE BY TRANSPOSING THE A MATRIX (AT) 
DO 61 I=1,M 
DO 61 J=1,N 

ea AT (J ,1) =A (I,J) 

C 

C USING STEPS (1), (2), AND (3) COMPUTE THE INVERSE 

e OP THE TRANSPOSE (AT) * WZIGHTED MATRIC (QLL) * 

C WETRIX A = OXX- 
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7) 


71 


8 1 


304 
B02 


B05 


306 
303 


507 


(7) COMPUTE AQ = AT * QOLL 


DO 71 I=1,N 

DO 71 J=1,M 

moun J) =Ol 

DO 71 K=1,M 

AQ (I,J) =AQ (I,J) +(AT(I,K)*QLL(K, J) ) 


(2) COMPUTE QXX = AQ * A 


DO 81 I=1,N 
DO 81 J=1,N 
SIS Ct pa ae 
DO 81 K=1,4M 
QXX (I,J) =OXX (I, J) +AQ (I,K) ¥A(K, J) 


(3) INVERT QXX MATRIX 


DO 307 K=1,N 

Bersd2 J=1,' 

TF (J-K) 304, 302, 304 

DeGiK yd) =O%X (K, J) /OXX (K, i) 
CONTINUE 

pec CMe) = 170 XX {K, *) 

DO 307 I=1,N 

IF (I-K) 305,307, 305 

DO 303 J=1,N 

IF (J-K) 306,303, 306 

Pend j= OX (L, J) —OXX (1, K) OX CUR, J) 
CONTINUE 
QXX(I,X)=-QXX (I,K) *OXX (K, K) 
CONTINUE 


i 





Cr Gy oO) 


Geren Cl aco 


io 1 


201 


301 


USING STEPS (4) AND (5), COMPUTE THE UNKNOWNS X 
BY MULT THE INVERSE QXX AND AQL. 


(4) COMPUTE AQL = AQ * EL 


DO 101 I=1,N 

AQL(I,1)=0. 

DO 101 K=1,M 

AOL(l, 1) =AOL (1, 1) +A (I, K) SEL (K, 1) 


(5) COMPUTE X = QXX * AQL 
DO 201 IT=1,N 


X (I, 1) =0. 
DO 201 K=1,N 
X (I, 1) =X(I,1) *QXX (I,K) *AQL (K, 1) 


(6) COMPUTE TEE RESIDUAL (V = A * X -EL) 


DO 301 I=1,M 

wae, 1) =0- 

DO 301 K=1,N 

V (I, 1)=V(I, 1) +A (I, K) *X (K, 1) 
DO 1 I=1,M 

ees 1) = vent) —EL (I, 1) 


COMPUTE THE STANDARD DEVIATION OF UNIT WEIGHT SIGMA 
Dre = NOMBER OF OBSERVATIONS 
DN - NUMBER OF KNKNOWNS 


SIGMA=0. 
DM=M 

DN=N 

BO 232 I[=1,m 


i 





YA AA 


OQ 


Ova a A 


382 SIGMA=SIGMA+YV (I, 1) **2 *QLL (1,1) 


446 


5 10 


SIGMA=DSQORT (SIGMA/(DM-DN) ) 


COMPUTE THE STANDARD DEVIATION OF THE ADJUSTED UNKNOWNS 
Tee Aas thE ELEMENTS CF THE COVARIANCE MATRIA. 


DO 446 I=1,N 
VAR(LI) =DSORT (QXX(I,1) *SIGMA**2) 


CONTINUE 
RETU RN 
END 


> a ee Pe ee eee ee eee eee eee ee eee ee ee ee ee eee ee eee eee ee ee eee ee eee ee ee ee 
SS ee eT ee ee ee ee EE ee 


Solves fOn CHE SHEMET MAJOR AND SEMI MINOR AXIS OF 
Pion OR BLLIPSE 


PMPLECET REALS (A-H,0-2) 

more tim woOk/7 A (30,30) , 50 (30, 1) ,OLL (30,30) ,AT(30, 30), 
remo sy, Ox X( 30,30) , AQE(30,10), X(30,1), V¥(30,1), 
eean(30), TITLE( 80) 


SUS = .5* (QXX(1,1) + QXX(2,2) + DSORT (DABS (QXX (1,1) - 
1 QXX (2,2) +4.0*0 XX (1,2) *OXX (2, 1)))) 
scr =O XX (1,1) * OXX(2,2) - DSORT (DABS (QXX(1,1)- 
1 QOXX (2,2) #4.0*9 XX (1,2) *OXX (2,1)))) 


SU = DSOQRT(SUS) 
SV = DSQRT (SVS) 
RETURN 


END 





Simul eNne GUESS DETERMINES THE BEST GUESS COORDINATES 
Bemoceuono EN SUBROUTINE DSORs. USE RIGHT SIDE RULE FOR 
SLALEON ORDER. 


IMPLICIT REAL*8 (A-H, 0-2) 
CCMMON /ISTAT/ STA (30,3),FXY(30) 


Dice eMtNE DESTANCS BETWEEN STATION 1 AND STATION 2 


Dower ((Sieame, 1) Sta (1, 1) ) **2+ (STA (2, 2) <STA (1, 2) ) **2) 


DETERMINE ANGLE ALPHA BETWEEN XO,YO/STA2/STA1 


ALPHA=DARCOS ( (STA(2,3) **2-STA (1, 3) **®2+D*®*2) / 


2s 0S tA (2, 37% D)) 


DEtonoeNGs ANGLE BROVO BETWEEN X-AXIS AND STA2-STAl 


BROVO=DARSIN ( (STA (1,2) -STA(2,2))/D) 


DETERMINE X AND Y LENGTH 


X=STA (2,3) *DCOS (ALPHA+BROVO) 
MSA (2,5) “DS IN (ALPHA +BROVO}) 


DETERMINE XO AND YO 


XO 


STA(2,1) +X 
STA (2,2) +Y 


< 
© 
it 
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RETURN 
END 


IMPLICIT REAL*8 (A-H,0O-Z) 

DATA E2,RKO,A /.00676 8658D0,. 99998D0,6378206.4D0/ 
DATA RKGEO,W1,W2 /0.048481368D0,0.11422D0,21.73607D0/ 
DATA W3,W4/5104.57338D0,6367399.689D0/ 

DATA RADSEC /.0000048481368111D0/ 


RADLAT = SECLAT * RADSEC 

SINLAT = DSIN(RADLAT) 

SIN2LA = SINLAT * SINLAT 

COSLAT = DCOS (RADLA?) 

COS2LA = COSLAT * COSLAT 

P = (CMER - SECLON) / 10000.0D0 
V = A / DSQRT(1.0D0 - E2*SIN2LA) 
TANCON = 1.0D0 - SIN2LA/COS2LA 

S = WY * (RADLAT - SINLAT*COSLAT/10.0D0**6 +* 
* (W3 -COS 2LA* (W2-W1*COS2LA) ) ) 
tiv= Sateen - CLAT 

ie) = RKGEOM*MGOSIAT * RKO IT 

73 = T2 * RKGEO / 2.0D0 

m= 73 * SINLAT 

BS= 13 M8 RKGEO * COS2LA / 3.0D0 
T6 = T5 * TANCON 

T7 = (4.0D0 + TANCON) *T5 *RKGEO * SINLAT / 4.0D0 
XCO = (T2 + (T6*P**2))*D + FEST 
YCO = (T7*¥P¥*uU) + (TU*D%K2) + 71 
RETURN 

END 


i 





> ee cee ce ee ee ee ee eee ee eee ee ee ee ee ee eee ee ee ee ee ee ee eee eee ee ee ee ee eee eee eee eee ee 
=_eo— oe oe ee ee ee ee ee ee ee ee ee eee eee es eee eee eee ee ee ee eee ee ee eee ee eee eee eee ee 


Pye el PR SALSS (A=H 70-2) 


DATA E2Z,A,SR /0-00676 86530D0,6378206.4D0, 
1 0.000004848136 8D0/ 
DATA W1,W2,W3 /0. 2468 2D0, 30.02335D0,5078.64977D0/ 


D = CLAT + YCO 

WO = 0. 15704998 1D0/10.0D0 **6 * D 
SINWO = DSIN(WO) 

COSWO = DCOS (WO) 

COS2WO = COSWO * COSHO 

PHI1 = WO + SINWO*COS YO/10.0D0**6 * 
* (W3+COS2 WO * (W249 1*COS2HO) ) 
PHI2 = PHI1 / 0.99998D0 

PHI3 = PHI2 / SR 


Q = (XCO - FEST) / 10.0D0**6 

V = A / DSORT(1.0D0 - E2*DSIN(PHI2) **2) 
T = DCOS(PHI2) * SR 

Cc = Vv * 0.99998D0 

T1 = 10.0D0**6 / (TC) 

T2 = (T1*10.0D0**6) / (2.0D0*C) 

T3 = (T2*10.0D0**6) / (3. 0D0*C) 


T4 = (T3*10.0D0**6) / (4. 0D0*C) 

DELLON = (T1-Q**2*T3* (2.0DO*DTAN (PHI2) **2+1.0D0) ) *0 
SECLAT = (((3-0DO*DTA N(PHI2) **2+5.0D0) *T4) *Q¥*2-T2) 
1 *Q**2*DSIN(PHI2) + PHI3 

SECLON = CMER - DELLON 

RETURN 

BND 
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LATITUDE 


(D=-M- 


Bio) 36 
36 36 
36 36 
ao 36 
36 36 


S) 


40.20 
42.41 
44.56 
46.84 
49.09 


DATA SET EXAMPLE 


Saee > POSLTION 


LONGITUDE 
(D-M-S) 


121 
120 
121 
en 
121 


52 
a2 
a2 
az 
52 


48.62 
48.16 
Gye 
46.98 
46.59 


V7 


PROGRAM GPBYLQ. 


OBSERVED LORAN RATES 
994 St Oa ye 


u=% 
(usec) 


275085) 9 
27508.81 
27508.83 
27508.76 
27 508.76 


(psec) 
42742.71 
42742.93 
42743.15 
§2743.26 
42743.48 





PROGRAM LORAN 


CY Cee) ee) 


ieee ikr REAL*S (A-H,O-Z) 
DIMENSION AXIS(13) ,RF (13) 


PRCGRAMN COMP UlES LINZ OF POSITIONS FOR LORAN-C USING 
THE SECONDARY FACTOR (SF) BASED ON SEAWATER EM MODEL. 


DIMENSION XSLAV1(1000),XSLAV2 (1000) 

DATA RHOSEC,PI, UNCOV, RN/2.062648062471D05, 
tous I26555090D0,299. 79 2458p0, 1.000338D0/ 
DATA XMEAN1,XMEAN2,VAR1,VAR2/9.00D0,0.00D0, 


1 0.00D0,00.0D0/ 


DATA XMEAN3,XMEAN4Y/0. 00D0,0.00D0/ 


C 


CCC CCC RK He A He He He He He ee ae he ee ee ote Heke eK OK OK ee kk 


DATA AXIS /6.3782064D06,6. 378388 


DOG P7Oes 47 SINS 5DL06, 


foe O1 5 D06,6. 3/78 16D 06,6.378249145D06,6.378165D06, 
yo 378160D06,6. 3/8165 D06,6.378145D006,6.3775634D06, 


peo. 5 76245D06,6. 3/8135 0D06/ 
DAeh RE 76. 35058 38 D06, 2.97 


DOCepe ato tacolzooDOZ, 


fee 202, 2.90247 16 / 427002, 2.93 465D02 ,2.9825D02, 


ee 02, 2.98 30D02,2.9825 D02, 6 


Beeeeoos DO2,2.98 26D02/ 
Cee ee* ELLIPSOID OPTION NUMBER 


1. CLARKE 1866 op 
C 2. INTERNATIONAL (HAYFORD) 9. 
C 3. BESSEL 1841 10. 
c GW. AND (AUSTRALIAN) ie 
c 5. 1967 REFERENCE es 


78 


oD O2 509006, 


Pane u Ry 
MARSHALL ISLAND 
NAVY 8D 

ATRY 

KRASSOWSKI 1940 






c 
c 
c 


6. CLARKE 1880 MOD bsiep WGS 1972 
7. SAO 


CCC CCC RK HK HM KK eK KK KK KKK 


e 
c 
c 


100 


Tr AO a ea Ga 


202 


201 


105 


Geta? = SLLISOSOID NUMBER (K) 
CC3-5 = NUMBER OF POINTS ALONG COAST (IREC) 


READ (5,100) K,IREC 
FORMAT(1X,12,1I4) 
TWOPI=2.*PI 

A=AXIS (K) 

F=1. /RF (K) 

IF(F.LT.3.D-3) F=(A-1./F) /A 


Peooeuaotek AND SLAVE STATIONS POSITIONS 
DHeerukS® RECORD IS THE NUMBER OF MASTER AND SLAVE 
STATIONS FOR THE PARTICULAR CHAIN. 


WRITE (6,202) 

FORMAT (181) 

WRITE (6,201) 

FORMAT (1H ) 

READ (5,105) ILATM,IMINM,RSECM,ILONM,ILMINM, RLSECM 
WRITE(6,105) ILATM,IMINM, RSECM, ILONM, ILMINM, RLSECH 
Fem rr (1%,03,1X,12,1£,FS. 2, 1X,14, 1X,12, 1X, FS.2) 
RMAST = ((IABS(ILATM) * 60 + IMINM) * 60 + RSECM) / 


Hens EC 


mentee bt. 0) RMAST = -RAAST 
Poodle = eA DS (LLONM) # 60 + LT LMINM) *60 +RLSECH) / 


1 RHOSEC 


Boe GONM GT. 0) RMASTL = TWOPE = RMASTL 


Wetton trot soLAVS STALLION AND CHANGE THe LATITUDES 
AND LONGITUDE INTO RADIANS. 


a9 





hennio, Wie Lia lS, (MINS, RSECS, [ LONS, ILMINS, RLSECS, 
1 DELAY! 
Soo, Very LOANS, LMINS, RSECS, TLONS, ILMINS, RLSECS, 
1 DELAY1 

101 Po eebneeeS VAP cd ky hoe 2, 1X, 4, 1X, 12, 1X, F5.2,F9.2) 
Romney i= ((2ABS (TLATS) = 60 + IMINS) ** 60 + RSECS) / 


1 RHOSEC 

IF(ILATS .LT. 0) RSLAV1 = -RSLAV1 

RSLAL1 = ((IABS (ILONS) * 60 + ILMINS) * 60 + RLSECS) / 
1 RHOSEC 


IF (ILONS .GT. 0) RSLAL1 = TWOPI - RSLAL1 


READ THE SECOND SLAVE STATION AND CHANGE THE LATITUDE 
AND LONGITUDE INTO RADIANS. 


QI AA YN 


READ (5,101) ILATS,ININS,RSECS,ILONS, ILMINS, RLSECS, 
1 DELAY2 

WRITE (6,101) ILATS,IMINS, RSECS, ILONS, ILMINS, RLSECS, 

1 DELAY2 

RSLAV2 = ((IABS (ILATS) * 60 + IMINS) * 60 + RSzZCS) / 

1 RHOSEC 

IF (ILATS .LT. 0) RSLAV2 = -RSLAV2 

RSLAL2 = ((IABS (ILONS) * 60 + ILMINS) * 60 + RLSECS) / 
1 RHOSEC 

IF(ILONS .GT. 0) RSLAL2 = TWOPI - RSLAL2 

WRITE (6,201) 


PT onoeaeecORDS POST TIONS AND OBSERVED KATES, THEN 
SOmr ures THn IOTAL DISTANCE. 


CY Oa A 


DOD 50 L=1ptREC SI 
ane omer) LLAT, LMIN ,RSEC,T LON, LEMIN, RLSEC, 
lsat et, RAT 32 


30 





mot  FORMAT(1X,13,13,F6.3,15,13,F6.3,3X%,2F10.2) 
Cc 
C CONVERT LAT AND LONG TO RADIANS. 


RPOST = ((IABS(ILAT) * 60 + IMIN) * 60 + RSEC) / 

1 RHOSEC 

IF(ILAT .LT. 0) RPOST = -RPOST 

RPOSTL = ((IABS(ILON) * 60 + ILMIN) * 60 + RLSEC) / 
1 RHOSEC 

IF(ILON .GT. 0) RPOSTL = TWOPI - RPOSTL 


SOrePULe DISTANCES AND AZIMNUTHS FROM THE OSSERVED POINT 


BAS rie 


QA AA A 


CALL INVER1(A,P,RPOST,RPOSTL, RMAST, RMASTL, FAZM, BAZM, 
1 DISTM) 

ODISTM = (RN * DISTM) / UNCOY 

CALL SECFAC (UDISTM,SF4) 


SLAV#1 


oun Ven h(a, >RPOST ,KRPOSTL,RSLAV1,RSLALI,FAZ1,5AzZ1, 
toed £> Tf 1) 

pee ot —= (RN * DEIST1) / UNCOV 

eVvote ser ac {UDEST1,SF 1) 


SLAVE2 
ServeenveR) (A, 2, RP Ost ,RPOSTL,RSLAV2Z, RSLALZ, FAZ2, BAZ2, 
ee iS T 2) 


UDIST2 = (RN * DIST2) / UNCOV 
Grat) SECFAC (UDIST?2,Sr 2) 
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COMPUTE THE RATES AND COMPARE [TO THE OBSERVED RATES 


© 


fiom UDTSt+ — UDISTN + SFI - SEM.% DELAY! 
DIFF1 = TDM1 - RATE 


moze Uuntol2 — UDISTM + SEZ — SPM + DELAY2 
Perr 2 = TDN2 - RATEZ 


COMPUTE THE LANE WIDTH IN METERS BASED ON EQUATION 
4.20 IN ELECTRONIC SURVEYING AND NAVIGATION - 
LAURILA, PAGE 94. 


AQAA A A 


BR1 = DABS(FAZ1 - FAZM) 
BR2 = DABS(FAZ2 - FAZM) 
WIDTH1 = (DIFF1 * UNCOV * 0.5) / DSIN(BR1 * 0.5) 
WIDTH2 = (DIPF2 * UNCOV * 0.5) / DSIN(BR2 * 0.5) 


| 
‘ 


Toten Oot TON OF VESSEL, COMPUTED RATE, OBSERVED 
RATE, AND THE DIFF BETWEEN THEM. 


CrCl el eo 


WRITE (6,200) ILAT,IMIN, RSEC,ILON, ILMIN, RLSEC,RATE1, 

DIFF 1,RATE2, DIF F2 
mgeFORMAT(1X,13,13,1X,F6.3,15,13,1X,F6.3,F12.2,F8.2, 

er i2 62, FS. 2) 

WRITE (6,201) 


c 
e PerveaN | AND XMEANZ ARE THE MSAN DIFFS BETWEEN THE 
C S@RBE URED RATS AND THE OBERSERVED. XKSLAV1I AND XSLAV2 
c Momecnn SHORED DIFFS. 
ic 
XMEANT = XMEANT + DIFF! 
XMEAN2 = XMEAN2 + DIFF2 
XMEAN3 = XMEAN3S + WIDTH1 


XMEANG = XMEANUY + WIDTH2 


OZ 





550 


600 


210 


211 


PmedO.3,° DIESZANCE IN METERS 


meio s,* DISTANCE IN METERS 


XSLAV1(I) = DIFF1 
XSLAV2(I) = DIFF2 
CONTINUE 


COMPUTE THE MEAN AND STANDARD DEVIATION 


XMEAN1T = XMEAN1 / IREC 
XMEAN2 = XMEAN2 / IREC 
XMEAN3 = XMEAN3 / IREC 
XMEANG4 = XMEANY Y/Y TREC 


DO 600 I = 1,IREC,1 
VAR1 = VAR1 + ((XSLAV1(I) - XMEAN1) **2) 
VAR2 = VAR2 + ((XSLAV2(I) - XMEAN2) **2) 


VAR1 = VAR1 / (IREC - 1.0) 
VAR2 = VAR2 / (IREC - 1.0) 
SD1 = DSQRT (VAR 1) 


SD2 = DSQRT (VAR 2) 


WRITE (6,201) 

WRITE (6,210) XMEAN1,SD1,XMEAN3 

WRITE (6,201) 

WRITE (6,211) XMEAN2,SD2,XMEANG 

WRITE (6,202) 

FORMAT(1X," SLAVE #1, MEAN = ',?10.3, 
‘ STANDARD DEVIATION = ', 


eS 
moma n( (X5" SLAVE #2, MEAN = ',F10.3, 
* STANDARD DEVIATION = ', 


Os 3) 
SLOP 
END 
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C 

c Peis eROULINS Wilh COM PULE THE SEA SECONDARY FACTOR 
e Wrpest = fOTAL DISTANCE 

C SF = SECONDARY FACTOR 

& 


IMPLICIT REAL*¥8 (A-H, 0-2) 


SCQOEPRPICIENTS 


QA A A 


123202323 
-~0.40758 
0.00064576813 


AQ 
A 
are 


2.741282 
-0.011402 
-000 32774815 


BO 
BY 
B2 


PED EST =<Gr. 537) SO TO 10 
ae B07 UTDEST) + B1 + ( B2 * UTDIST) 
GO TO 20 

10 SF = ( AO / UTDIST) + Al + ( A2 * UTDIST) 


20 SONTINUE 
RETURN 
END 


=e ee eee ces ee eee eee ee ee ee ee ee eee ee eee ee ees ee ees ee es ee eee ee eee ae 
_—en oe ee ee ee ee eee ese eee eee ee eee es ee eee ee eee es ee ee eee eee eee ee 


Sue OUGLENE ENVSRI (A, INV, GLAT?1,GLON1,GLAT2Z,GLON2,FAZ, 
1 BAZ,S) 


8u 





CY OY OY OY OY OY Ora Ye Ol) Cle era 


ete SOLUTZON OF THE GEODETIC INVERSE PROBLEM AFTER 

are To VENGENTY MODIFIED RAINS FORD'S METHOD WITH HELMERT'S 
eee PLELPTICaL TERNS. EFFECTIVE iN ANY AZIMUTH AND AT 
weer any DESTANCE SHORT OF ANTIPODAL STANDPOINT/FOREPOINT 
*x**% MUST NOT BE THe GEOGRAPHIC POLE 


Sitmr (> THE SEMI-MAJOR AXIS OF THE REFERENCE ELLIPSOID 
See (NY 15S THE FLATTENING (NOT RECIPROCAL) OF THE 

SSF REFERNECE ELLIPSOID LATITUDES AND LONGITUDES IN 

¥*¥* RADIANS POSITIVE NORTH AND EAST FORWARD AZIMUTHS AT 
mee BOTH POINTS RETURNED IN RADIANS FROM NORTH GEODESIC 
oer SOP OLANCE S RETURNED IN UNITS OF SEMI-MAJOR AXIS A 


*** DROGRAMMED FOR CDC~-6600 BY LCDR L.PFEIFER NGS 
*x** ROCKVILLS MD 18FEB75. NOl meno FOR LSM SYSTEM 360 
**x BY JOHN G GERGEN NGS ROCKVILLE MD 7507. 


IMPLICIT REAL*8 (A-H, 0-2) 
DATA =PS/0.5D-13/, PI/3. 1415926535898D0/ 
TWOPI=2.*PI 
R=1.-FINV 
TUI=R*DSIN(GLAT 1) /DCOS(GLAT1) 
TU2=R¥*DSIN(GLAT 2) /DCOS (GLAT2) 
CU1=1./DSQRT (TU 1*TU1+ 1.) 
SU1=CU1*TU1 
CU2Z=1./DSQRT (TU 2*TU2+ 1.) 
S=CU 1¥*C U2 
BAZ=S*TU2 
PAZ=BAZ*TU1 
X=GLON2-GLON1 

100 SX=DSIN (X) 
CX=DCOS {X) 
7 1=CU2*S X 
TU 2=BAZ-SU1*CU2 *CX 
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SY=DSQRT (TUI*TU 14+TU2* TU2) 
CY=S*CX+FAZ 
Y=DATAN2(SY,CY) 
SA=S*SX/SY 
C2A=-SA*SA+1. 
CZ=FAZ+FAZ 
EF (C2A.GT.0.) CZ=-CZ/C 2A4CY 
B=CZ*CZ*2.-1. 
C= ((-3. *C2A44 .) *FINV+ 4.) *C2A*¥FINV/16. 
p=X 
X= ((E*CY*C+CZ) *SY*C+Y) *SA 
X= (1.-C) *X*FINV4+GLON2 -GLON1 
IF (DABS (D-X).GT.EPS) GO TO 100 
FAZ=DATAN2 (TU1, TU2) 
IF (FAZ.GE.TWOPI) FAZ=FAZ-TWOPI 
IF(FAZ.LT.0.D0) FAZ=FAZ+TWOPT 
BAZ=DATAN 2(CU1*SX, BAZ *CX-SU1*CU 2) +PI 
IF(BAZ.GE.TWOPI) BAZ=BAZ-TWOPI 
IF (BAZ.LT.0.D0) BAZ=BAZ+TWOPI 
X=DSQORT ((1./R/R-1.) #C 2A+1.) #1. 
Mori 2. 1/7 x 
eae x 
C= (X*xX/4.4+1.)/C 
D=(0.375*X*X- 1.) *X 
X=E*CY 
S=1.-E-E 
S=(( ( (SY*SY*4.- 3.) *S*CZ*D/6.-X) ¥D/U.+CZ) *SL*D+EY) ¥CKAER 
RETURN 
END 

SENTRY 
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DATA SET EXAMPLE = PROGRAM LORAN 


Spoiler 5) POSETION OpsSe hue) hORAN ORAT SS & ERRORS e) 
LATITUDE LONGITUDE S94U0-Y oY oo40-W Ee 
Co =4'— S) (D=iN=5) (psec) (usec) (asec) (psec) 


36 43 45.800 121 55 27.160 42789.34 -0.49 16294.04 ~-1.96 
some 632400 121 55 32.380 42791.13 -0.38 - 16293.46 -1.10 
PomeeneioOm lara 37.390 42793204 -0.38 16292.73 -0.99 
memes 76490 127 55 46.950 42795.13 -0.58 16292.03 -1.96 
Somes 53.260 121 55 57.719 42796.93 -0.51 16291.43 -1.27 
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PROGRAM LOPLC 


PeeGnaAn COMPUTES LINE OF POSITIONS FOR LORAN-C USING 
PoioeoSCONDARY FACTOR (SF) AND ADDITICNAL SECONDARY 
ByCtOnmmtnasr) THE SP AND ASF ARE BASED ON 

Semen PNPiRECAL TD GRID. 


CY OY A OY Cte oO 


IMPLICIT REAL*8 (A-H, 0-2) 
DIMENSION PHI (100) ,ELON (100) ,DISTM (100) ,DIST1(100), 
1 DIST2(100) , AXIS (13) , RF(13) ,AZ1 (100) ,AZ2(100) , AZM (100) 
DATA RHOSEC, PI, UNCOV/ 2.06 2648062471D05, 
1 3.1415926535898D0 ,299.693D0/ 
Cc 
CCT CCC RK RK RM KR KK KK KK KKK KKK KK K 
DATA AXIS/6.3782064D06,6. 378388 D06,6.377397155D06, 
1 6.37816 D06,6. 37816D 06,6. 378249145D06,6.378165D06, 
2 6.378166D06,6. 378165 D06,6.378145D06,6. 3775634006, 
3 6.378245D06,6. 378135 0D06/ 
DATA RF/6. 35658 38 D06,2.97 D02,2.9915281285D02, 
1 2.9825D02, 2. 98247 167427D02,2.93465D02,2.9825D02, 
2 2.983 D02,2.983D02,2.9825 D002, 6.3562569D06, 
3 2.983 D02,2.98 26D02/ 
C*****ELLIPSOID OPTION NUMBER 


c 1. CLARKE 1866 SeeeoneU RY 

Cc 2. INTERNATIONAL (HAY FORD) 9. MARSHALL ISLAND 
c 3. BESSEL 1841 10. NAVY 8D 

c G4. AND (AUSTRALIAN) Pema Re, 

c Seaeioo? REFERENCE 12. KRASSOWSKI 1940 
© 6. CLARKE 1880 MOD \e2rnGsS. 1972 
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e 
e 


io AO 


CCC CC CO I IK HK eK KK KKK eK OK 


. 
c 
3 


100 


105 


Oo Cy er 


101 


CC1i-2 = ELLIPSOID NUMBER (K) 
CC3-5 = NUMBER OF POINTS ALONG COAST (IREC) 


READ (5,100) K,IREC,ISTA1, 1STA2 
BORMAT (1X, 12, 04 , 223) 
TWOPI=2.¥*PI 

A=AXIS (K) 

F=1. /RF (K) 

IF (F.LT.3.D-3) F=(A-1./F)/A 


ero MASTER STATEON, POSITIONS. 
eens) LLATM, EMI NM,RSECM, LLONM, ILMINN, RLSECN 


Opa 1k 23, 4%,02,1% ,F5.2,1X,i8,1X%,12, 1X, F5-2) 
RMAST = ((IABS(ILATM) * 60 + IMINM) * 60 + RSECM) / 


ten HOS EC 


IF(ILATM .LT. 0) RMAST = -RMAST 

RMASTL = ((IABS (ILONM) * 60 + ILMINM) *60 +RLSECM) / 
RHO SEC 

IF(ILONM .GT. 0) RMASTL = TWOPI - RMASTL 


Poni Pinot SLAVS STAPION AND CHANGE THE LATITUDE 
Poo NGONGIIUDE INTO RADIANS.DELAY iS THE CODING DELAY 


Migeotac to tHe OFFSET IN MICRPOSECONDS. 


Peer VO 1) LLATS, 2Mi NS, ,RSEGS, LTLONS, ILMINS, RLSECS, 


feveLbAyY|, SLAS% 


emer So nlIk Gl 2, tk, ho. 2, Xe l 4, 1X,12,1X,F5.2, 


eno. 2,P7. 3) 


RSLAV1 = ((IABS (ILATS) * 60 + IMINS) * 60 + RSECS) / 


1 REOSEC 
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IF(ILATS .LT. 0) RSLAV1 = -RSLAV1 

RSLAL1 = ((IABS (ILONS) * 60 + ILMINS) * 60 + RLSECS) / 
1 RHOS=C 

IF(ILONS .GT. 0) RSLAL1 = TWOPT - RSLAL1 


HoeAvelis SSCOND SLAVa STATION AND CHANGE THE LATITUDE 
AND LONGITUDE INTO RADIANS. 


Cry Or) 


READ (5,101) ILATS, IMINS,RSECS,ILONS,ILMINS,RLSECS, 

1 DELAY2,BIAS2 

RSLAV2 = ((IABS(ILATS) * 60 + IMINS) * 60 + RSECS) / 

1 RHOSEC 

IF(ILATS .LT. 0) RSLAV2 = -RSLAV2 

RSLAL2 = ((IABS (ILONS) * 60 + ILMINS) * 60 + RLSECS) / 
1 RHOSEC 

IF(ILONS .GT. 0) RSLAL2 = TWOPI ~ RSLAL2 


READ COAST POINT LAT AND LONG AND CONVERT TO RADIANS. 
STORE LAT IN PHI AND LONG IN ELON. 


Aa aA NM 


DO 500 I=1,IRzC,?1 

READ (5,102) INUM,ILAT,IMIN, RSEC,ILON,ILMIN, RLSEC 
fem FORMAT(1X,14,24,13,F7-.3,15,13,F7-3) 

PHI(I) = ((IABS (ILAT) * 60 + IMIN) * 60 + RSEC) / 


leno Sec 

me(@LAT .<LiT. 9) PHI(I) = -PHI (7) 

peer) —s( (LABS (LLIN) * 60 + TLMEN) * 60 + RLSEC) / 
lero S sc 


IF(ILON .GT. 0) ELON(I) = TWOPI - ELON(TZ) 


Semeur e DISTANCE FROM COAST POENT TO MASTER AWD 
Save StATRONS. SET UP COMPUTATION AND CALL INV 


Ore a 


P1 = PHI(T) 





£1 = ELON (I) 

CALL INVER1(A,F,P1,£1,RMAST,RMASTL,AZF,AZB,S) 
DISTM(I) = S 

AZM(I) = AZF 

CALL INVER1(A,F,P1,E1,RSLAV1,RSLAL1, AZF,AZB,S) 
iment) = Ss 

AZ1(I) = AZP 

CALL INVER1(A,F,P1,£1,RSLAV2,RSLAL2, AZF,AZB,S) 
plete (i) = s 

Maat) = AZIE 

WRITE (6,210) I, INUM,I LAT, IMIN,RSEC,ILON,ILMIN, RLSEC, 
4 AZM(I) ,AZ1(I) , AZ2 (I) 

mORMAT( 1% ,24,04 ,14 13,77. 3,15,13,F7.3,3F10. 3) 
500 CONTINUE 


bo 
—- 
© 


Paves NUMBER OF RECORDS AND THEIR POSITIONS. 

Sere ut eens EBO@PAL DISTANCE OVER LAND AND SEA. 
Pio ieee nOnAL DEISTANGE FROM BASTER TO DATS SCOINT. 
Peo bona. DISTANCE FROM SLAVE TO DATA .POINT. 


COMPUTE THE FORWARD AND BACK AZIMUTHS 
faye AND FPAZS = FORWARD AZINOTH TO THE MASTER AND 
Save ofA LON. 
ERoGueANOeSASS = BACK AZEMUTH TO THE MASTER AND 
Soave ATION. 


CG Gr) Or tee Cl kee ad 


READ (5,103) JREC 
FORMAT (1X,24) 


Q —« 
© 
LJ 


BONo50 L=4,0REC,7 
READ (5,104) LLAT ,LUIN, RSEC,ILON, ILMIN, RLSEC, RATE1, RATED 
104 meemar(1X,13,13,F6.3, 15,13,F6.3,F9.3,F9.3) 


. CONVERT LAT AND LONG TO RADIANS. 
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Cy ey el er 


RPOST = ((IABS(ILAT) * 60 + IMIN) * 60 + RSEC) / 

1 RHOSEC 

IF({ILAT .LT. 0) RPOST = -RPOST 

RPOSTL = ((IABS(ILON) * 60 + ILMIN) * 60 + RLSEC) / 
1 RHOSEC 

IF(ILON .GT. 0) RPOSTL = TWOPI - RPOSTL 


COMPUTE DISTANCES AND AZIMUTHS FROM THE OBSERVED 
Peryr TO Tht DLNTERPOLATED SHORE POINT AND WRITE. 


isch =. 7 

awe otOone nr (ERE C,RrOS TVRPOSIL, RMNAST,RMASTL,UTIDISA, 
Pauses tum ES l, FAZM,B AZM 

Sxeiwoece nc (UTDISt,0S DIST,ULDIST, BAZM,ISTA, SFM} 


Grime ouORPT (TREC, XPOST,RPOSTL,RSLAV1I,RSLAL1,UTDIS1, 
imooots. ,J ELDEST, FAZ1,8 AZ 1) 

Prblbecheeae (UTD 151,USDIST,ULDIST, BAZ1,ISTA1,SF1) 
SweemmonenrPre ent C,RPOS T,RPOSTL,RSLAV2,RSLAL2Z, UTDIS2, 
iu oD EST, ULDIST, FAZ 2,3 AZ 2) 

Seo Cr ACU Dp IS2,USDiST ULDIST, BAZZ,1ISTA2,SF2) 


Sere Urs GHE RATES AND COMPARE TO THE OSSERVED RATES 


ieee prst - ULDISM + SFI - SFM + DELAY1 + BIAS? 
Pieper = TDMT - RATE 


eee bis2 = ULDESMet SF2Z2 - SFM + DELAY2Z + BIAS2 
mune = EDMZ = RATEZ 


BOmeurte THE EANE WEDTH IN METERS BASED ON EQUATION 
PecomrN eG LEGLRONIC SURVEYING AND NAVIGATION - LAURILA, 


we 
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BR1 = DABS(FAZ1 - FAZM) 

Mater iosGre PI) BR1 = TWOPT = BR1 

BR2 = DABS(FAZ2 - FAZM) 

IF(BR2 .GT. PI) BR2 = TWOPI - BR2 

WIDTH1 = (DIFF1 * UNCOV * 0.5) / DSIN(BR1 * 0.5) 
WIDTH2 = (DIPF2 * UNCOV * 0.5) / DSIN(BR2 * 0.5) 


allo THeeePOSITION OF VESSEL, OBSERVED RATES, 
AO LHe Dinr ERENCES BALWEEN THEA. 


Yr S| 


eeeenO, 20008 LoAl, Imi NO hSke, Tt LON, [LDNIN,PLSEC,RAT=1, 
mem fi RATEZ, DIFF2 
200 ree, i 2, oe 65, 24, 2S, br). 3, F10.2,F6.2,F 10.2, F6. 2) 


C 
@ RieaNTeAND XMEANZ ARE THe MEAN DIFFS BETWEEN THE 
c COMPUTED RATE AND THE OBERSERVED. XSLAV! AND XSLAV2 
C Pace bas SLORme Dir kS. 
c 
XMEANT = XMEANT + DIFP! 
Mane = XMBANZ + DIP r2 
XMEAN3 = XMEAN3 + WIDTH! 
AMEANG = XMEANGY + WIDTH2 


XSLAV1(I) = DIFF 
XSLAV2(I) = DIFF2 
550 CONTINUE 


C 
E COMPUTE THE MEAN AND STANDARD DEVIATION 
C 

XMEAN1 = XMBAN1 / JREC 

XMEAN2 = XMEAN2 / JREC 

XMEAN3 = XMEAN3 / JREC 

XMEANG = XMEAN4S / JREC 





HOmegO r= 1, JREC,1 


VAR1 = VAR1 + ((XSLAV1(I) - XMEAN1) **2) 
600 VAR2 = VAR2 + ((XSLAV2(I) - XMEAN2) **2) 
C 

VAR1 = VAR1 / (JREC - 1.0) 

VAR2 = VAR2 / (JREC - 1.0) 
C 

SD1 = DSQRT (VAR 1) 

SD2 = DSQRT(VAR2) 
C 


WRITE (6,201) 
WRITE(6,210) XMEAN1,SD1,XMEAN3 
WRITE (6,201) 
WRITE (6,211) XMEAN2,SD2,XMEAN4 
WRITE (6,202) 
210 FORMAT(1X,' SLAVE #1, MEAN = ',F10.3, 


[le SLANDARD DEVIATION = °, 

Peeetue so, DESTANCE [TN METERS = ',F10.3) 
211 memeam gia, OLAVE 2, MEAN = ', F10.3, 

1 ' STANDARD DEVIATION = ', 

Peo. Spee STANCE IN METERS = * 7510.3) 


WRITE (6,201) 
201 PORMAT (1H ) 
202 FORMAT(1H1) 


SOP 
END 
Cee 22 s2s22225 225 222225 225255 S255 2225555 e555 5555555 52255 52=ee=== 
See UlONG SaOnRehGrns C,RPOST,RPOSTL, SG90UT, RCONTL,UDIST, 
1 USD,ULD, FA, 3BA) 
ee i ee a Se ee ee ee 
c Seewouriwe SHORPI WELL SELaCT A POINT ALONG THE SdAORE 
c memes so OUTLINED £xO8 NORTH BY SELECTED POINTS wWiTH 






NAAN AN 


Crtr YY OF A Oo Te! 


KNOWN LATITUDES AND LONGITUDES. THE SHORE POINT IS 
INTERPOLATED BETWEEN TWO KNOWN POINTS USING THE TOTAL 
DESTANCE BETWEEN THE POSITION AND THE CONTROL STATION 
Mimelis AGEOUIH BETWEEN THE SHORE POINTS AND THE 
RECEEVERS POSITION. 


IMPLICIT REAL*8 (A-H, O-Z) 

DIMENSION FAZM(100) 
COMMON/SHORE/PHI(100) , ELON (100) , UNCOV,RN,A,F 
DATA PI/3.14159 265358 98D0/ 


TWOPT = 2.0 $e PI 


CALL INVERT (A,F,RPOST ,RPOSTL,RCOONT, RCONTL, FA, BA, RDISTT) 
ODIST = (RN * RDISTT) / UNCOV 


RCOMP = 99999.99 

DO 10 J=1,IREC, 1 

Pae= PBT (J) 

Bl = £LON (J) 

CALL INVER1(A,F,P1,£1,RCONT, 3CONTL, AZF,AZB, RDISTL) 
CALL INVER1(A,F,RPOST,RPOSTL, P1,£1, FAZ,BAZ, RDISTS) 
USD = (RN * RDISTS) / UNCOV 

ULD = (RN * RDISTL) / UNCOYV 

FAZM (J) AZF 


Cte DIF PERENCE BEP@EEN THE TOTAL DISTANCE 

eves! jean) CHE SUMMATION OF THE DISTANCE OVER THE WATER 
CoD) AND THE DISTANCE OVER THE LAND (ULD). IT THE 

Pe eANC Gels LESo GTEAN RCOMP, UPDAT RCOMP AND JSTAa. 

tee aeto THE CLOSESI POINT ALONG THE SHORELINE WHICH I5 
NEAR THE 2M PROPAGATION PATH. 





QaAaAAN = 


EDyee = DABS(UDZST - (USD +-ULD)) 
Pro urbe Gls —RKCOMP) GO TO 10 
REGnNe = FDIFF 

JSTA = J 

CONTINUE 


NOW DETERMINE THE INTERPOLATED LATITUDE AND LONGITUDE 
[oem “OL Nt USING AZT MUTE PERCENTAGE. 


IUPPER = JSTA - 1 

ILOWER = JSTA + 1 

AZMU = FAZM(IUPPER) 

AZML FAZM (I LOWER) 

IF (AZMU .GT. FA .AND. FA .GE. FAZM(JSTA)) ICH = IUPPER 
LP(FAZN(JSTA) .GE. FA .AND. FA .GT. AZML) ICH = ITLOWER 


RADJ = 1.00 - DABS((FAZM(ICH) - FA) / (FAZM(ICH) - 


1 FAZM(JSTA) )) 


RNWPHI = PHI(JSTA) + ((PHI(ICH) - PHI{(JSTA)) * RADJ) 
IF (ELON (JSTA) .LE. ELON(ICH)) RADJ = 1.0 - RADJ 
RNWELN = ZBLON(SSTA) + ((ELON(ICH) - ELON(JSTA)) * RADJ) 


CALL TODMS(RNWPHIT,IDG ,AIN, SEC) 
RHOLD = TWOPI - RNWELN 
CALL TODMS (RHOLD,IDGL,MINL,SECL) 


CALL INVER1(A,F,R3NWPHI, RNWELN, RCONT, RCONTL,AZ,BZ,RDISTL) 
CALL INVER1(A,F,RPOST,RPOSTL,RNWPHI, RNWELN,AZ,5Z, RDISTS) 
USD = (RN * RDISTS) / UNCOY 

ULD = (RN * RDISTL) / UNCOY 

RETU RN 

END 


ee ees ces es ee ee ee ees ee es es es ee eee ee ee ee ee ee ee ee ee ee eee ee eee ee ee es ee eee ee ee ee ee ee ee ee ee 


96 





OO AY Ga a a 1 a a ae ae 


SeneULIink SECFAC(UIDIST, USDIST, CLDIST,AZL,ISTA,SF) 


_—anP =P ee eee ee eee eee Eee eee ee eee eee ee ee Se See Se ee ee ee eee ewe ee eee ee ee ee See ee eee ee ee ee ee Se ee eee ee ee Se oe 


THIS ROUTINE WILL COMPUTE THE LAND/SEA SECONDARY FACTOR 
UrDist = TOTAL DISTANCE 
USDIST = DISTANCE OVER THE SEA WATER PATH 
ULDIST = DISTANCE OVER THE LAND PATH 
AZI = AZIMUTH FROM NORTH. 
MASTER = 1 (ISTA) 
W = 2 
X = 3 
ie = 4 


IMPLICIT REAL*8 (A-H, 0-2) 

S1 = (.795 / USDIST) + 0.439 + (.00245 * USDIST) 
IF(UTDIST .GT. 580) GO TO 10 

S2 = (3.188 / UTDIST) - 0.594 + (.000329 * UTDIST) 


GO TO 20 
10 S2 = (128.8 / UTDIST) + 0.187 + (.000652 * UTDIST) 
20 IF(USDIST .GT. 540) GO TO 30 
$3 = (3.188 / USDIST) - 0.594 + (.000329 * USDIST) 
GO TO 40 
30 See it26.e 7 USdIsr) + 9.187 + (.090652 « uSsDIsT 
40 IF (ULDIST .GT. 540) GO TO 50 
S4 = (3.188 / ULDIST) - 0.594 + (.000329 * ULDIST) 
GO TO 60 
50 Si = (128.8 / ULDISE) + 0.187 + (.000652 * SJLDIST) 
60 $5 = 1.428 + (.00158 * UTDIST) 


S6 = 1.428 + (.00158 * ULDIST) 
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TAZI = 2.0 * AZI 
Comino, (7/0,80,90,100) , ISTA 


70 RHOLD = (1. 010*DSIN (AZI)) - (.196*#DCOS (AZT) ) 
1 - (.893*DSIN(TAZI)) - (.355*DCOS (TAZI) ) 
GO TO 200 

80 RHOLD = (.323*DCOS (AZI)) - (.711*DSIN (TAZI)) 
GO TO 200 

90 RHOLD = (942 *DCOS (TAZI) ) 
GO TO 200 

100 RHOLD = (.588*DSIN(TAZI) ) 

Cc 

200 S5 = S5 + RHOLD 


S6 = S6 + RHOLD 
Soares = oS £ S6 = St + S2 + S353 - Sh) 
mer 46,500) S1,52,53 ,S4,94,S5,S6,SF 

500 FORMAT (1X,7F15. 5) 


RETURN 
ED 
Cas 2222222 222 55555555 SS2525 2255525555 555555 S55 S55 5 5555552552 
SueROULINE ENVERI{(A, F INV, GLAT1,GLON1,GLAT2,GLON2,FAZ, 
ne BAZ »S) 
Ore a a a a ee eee 


See Appendix B for subroutine INVER. 
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meorac DOlntS around Monterey Bay, California. 
are used to interpolate geodetic point 


of the distance over 


NO 
- 
2 
5 
4 
> 
6 
7 
8 


3 
10 
11 
12 
eS 
14 
is 
16 
A? 
18 
cps 
20 
Zt 
ae 
Z3 


BIG 
BiG 
36 
36 


2 
~ 


36 
Z:0 
36 
Sie 
210 
36 
36 
36 
36 
36 
a5 
36 
36 
36 
36 
36 
36 
36 


LALTEIUDE 


oy 
70 
Sih 
a 
57 
1 
58 
oo 
56 
5 
3 
49 
47 
46 
44 
41 
31S) 
a 
56 
36 
56 
a7 
aC 


18.606 
59.204 
05.076 
49.538 
17.949 
08.589 
32.140 
01.498 
46.115 
38.140 
13.806 
38.384 
39.241 
Dic aou 
DOT 
145439 
17.211 
Si 125 
23.446 
C020 
24.782 
licaao 4 
00.300 


For 


land and the distance 


LONGITUDE 


122 
122 
122 
| 
121 
121 
121 
121 
24 
AZ 
121 
127 
121 
121 
121 
121 
121 
et 
121 
en 
121 
V2 
24 


05 
03 
01 
01 
58 
57 
55 
53 
52 
51 
49 
U7 
47 
47 
47 
48 
49 
50 
51 
52 
53 
54 
55 


B05 Zo 
01.8177 
31.701 
C7557 
iga5 350 
07.288 
10.083 
572390 
Zee 
2453599 
46.743 
48.3895 
10.818 
ao. 
D2 a4 1c 


32.602 


ale eis ee 
a liawees 
34.833 
[705579 
48.453 
171.628 
Is 550 


tHe pOLns 
computation 


OVEE Sea. 





BATA SET EXAMPLE - PROGRAM LOPLC 
Shape’ Ss POSLTION 


36 
36 
36 
36 
36 


LATITUDS 
(D=ti= S) 


43 45.800 
G4 3.400 
44 21.180 
44 37.490 
44 53.260 


LONGI 
(D=t=S) 


121 
121 
121 
AY 
121 


25 
55 
25 
22 
2 


TUDE 


27. 160 
32. 340 
SoG 
46.950 
Soe sO 


(psec) 
42789. 34 
42791.13 
42793.04 
42795. 13 
42796.93 


100 


(NSEC) 
1.174 
2S 
lec 
1.05 
Nod 


(NSEC) 
16294 .04 
16293 .46 
GZS 
loZzo 2.03 
16291.43 


OBSERVED LORAN RATES & ERRORS Ne) 
a5 0 BY 9940-W mY 


(us ec) 
re 32 
ieee 
1. 40 
t53 
1.11 





OQ OA AGF Oye 


C 


PROGRAM LORTAB 


PROGRAM GCOMPUTES LINE OF POSITEON FOR LORAN-C USING 
SP SALT WATER CORRECTION FACTOR AND DMAHTC CALCULATED 
Cher LeiD OBSERVED ASF CORRECTIONS. 


IMPLICIT REAL*8 (A-H, 0-2) 

DIMENSION AXIS(13) ,RF (13) 

DIMENSION XSLAV1(1000) ,XSLAV2 (1000) 

DATA RHOSEC, PI, UNCOV, RN/2.062648062471D05, 
1 3.1415926535898D0,29 9. 79 2458D0 ,1.000338D0/ 
DATA XMEAN1,XMEAN2,VAR1,VAR2/0.00D0,0.00D0, 
1 6.00D0,00.0D0/ 

DATA XMEAN3,XMEAN4/O.00D0,0.00D0/ 


CCC CC CR RRR RK KKK HK Re He KK KK KK KKK KK 


DATA AXIS/6.3782064D06,6. 378388 POG 62519597 oop 0o, 
ion oor DUG, oO. 2/6 16D06,6.378249145D006,6.3738165D06, 
mmo 6 1 GON0G, 6.379165 D0G,6.378145D06,5.3775634D06, 
pmo. 76245D06,6. 378 135 0D06 / 

Bia, RE /6 63560838 DO6, 2.97 MiGs eo 2OoNe Ss oDO2, 
Meee 2500 2,2. 99247 16 7 427D02,2.93405D02,2.9825D02, 
peo 2,2. 76 3D02,2.9825 DOZ, 6. 3562569D06, 

Emcee OS DO2Z, 2-98 260027 


CeeetesLLIPSOID OPTION NUMBER 


C 
c 
e 


1. CLARKE 1866 Ss MERCURY 
2. INTERNATICNAL (HAY FORD) Sen anon wht) LSLAND 
Be SESSEL 1841 LOS eerANY 8D 


101 





Cr Cr eGY er ©) 


CCC CCC eK EK ee He He eC oe eR A KE eK EK KO OK OK 


c 
c 
C 


100 


QA AA NAN 


ZO2 


201 


105 


4. AND (AUSTRALIAN) 
5. 1967 REFERENCE 
6. CLARKE 1880 MOD 


Te SAO 


i 
Va. 
toc 


Geq—=2 = ELLIPSOID NUMBER (K) 
GCC3=5 = NUMBER OF POLENTS ALONG COAST (TREC) 


READ (5, 100) 


FORMAT (1X,I2,1I4,213) 
TWOPI=2.*PI 


A=AXIS (K) 


F=1. /8F (K) 


IF (F.LT.3.D-3) 


Keene yeeD ily m2 


F=(A-1./7) /A 


Ray 
KRASSOWSKI 1940 
WGS 1972 


Pav ewnnost anh AND SLAVE STATIONS POSITIONS 
Tiger atol RECORD 15 THE NUMBER OF MASTER AND SLAVE 


SeaelONs FOR THE PArT ICULAR CHAIw. 


WRITE (6,202) 
FORMAT (1H1) 
WRITE (6,201). 
FORMAT(1H ) 


READ (5, 105) 


WRITE (6,105) 


RMAST = 
1 RHOSEC 
IF (ILATM 
RMASTL = 
1 REOSEC 
IF (ILONM 


( (TABS (ILATM) 


m 
ewe 


({IABS (ILONN) 


BAG ig 


eA ee eh oe ON MN, LLMINA, RLS ECM 


Prien kote, SLONt, LLNINM, RLSECM 
ere, opeiN yg l2,1k,l os 25 1k, 4, 1X,12, 1%,75.2) 


0) RMAST = 


0) RMASTL = 


102 


Rao. 


TWOPT 


‘OOo etn) = 60 + RSECM) / 


TecUeoe or waNiyeto0 +ALSECM) 7 


RMASTL 





c Poe nimer Rot SLAVE STARTOH AND CHANGE THE LATITUDE 
c AND LONGITUDE INTO RADIANS. 


READ (5,101) ILATS, IMINS,RSECS,ILONS, ILMINS, RLSECS, 
1 DELAY’ 
WRITE(6,101) ILATS,IMINS,RSECS, ILONS,ILMINS, RLSECS, 
1 DELAY1 

ame FORMAT (1X%,13,1X,12,1X ,F5. 2, 1X,I4, 1%,12, 1%, F5.2, 79.2} 
RSLAV1 = ((IABS(ILATS) * 60 + IMINS) * 60 + RSECS) / 
1 RHOSEC 
MehATS .bT. 0} RSLAVY1 = -RSLAV1 
RSLAL1 = (({IABS(ILONS) * 60 + ILMINS) * 60 + RLSECS) / 
1 RHOSEC 
IF(ILONS .GT. 0) RSLAL1 = TWOPI - RSLAL1 


Pea era be s2eOND SLAYE STATION AND CHANGE THE LATITUDE 
AND LONGITUDE INTO RADIANS. 


CV OY OYA 


READ (5,101) ILATS,IMINS,2SECS,ILONS,iILMINS, RLSECS, 

1 DELAY2 

WRITE (6,101) ILATS,IMINS,RSECS,ILONS, ILMINS, RLSECS, 
1 DELAY2 

RSLAV2 = (({IABS(ILATS) * 60 + IMINS) * 60 + RSECS) / 
1 RHOSEC 

IF(ILATS .LT. 0) RSLAV2 = -RSLAV2 

RSLAL2 = ({TABS(ILONS) * 60 + ILMINS) * 60 + RLSECS) / 
1 RHOSEC 

memmEONS GT. 0) RSLALZ = THOPT - RSLAL2 

WRITE (6,201) 


Peaoe ine RECORDS POSITIONS AND OBSERVED RATES, THEN 
memrure THE TOTAL DISTANCE. 


Cr OG) 


DO 550 I=1,IRzC,1 
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Rowe >, (04) TLAT SEMEN ,RSeC, LLON, [LLMIN,RLSEC, 
1 RATE, RATE2 
104 ROMS, IS ,6O.s, L5,L3,6O.5, 5X, 2F 10.2) 


GS 


C CONVERT LAT AND LONG TO RADIANS. 
c 
RPOST = ((IABS(ILAT) * 60 + IMIN) * 60 + RSEC) / 
1 RHOSEC 
IF(ILAT .LT. 0) RPOST = -RPOST 
RPOSTL = ((IABS (ILON) * 60 + ILMIN) * 60 + RLSEC) / 
1 RHOSEC 


PMeerbON 2Gi. 0) @REOSTL = TWOPIT - RPOSTL 


Siicu re PESPANCES AND AZENOTHS FROM THE OBSERVED POINT 


MASTER 


CY I "Cy Ol el 


Crete ENVER?(A,ZP,RPOST ,RPOSTIL, RMAST,RMASTL, FAZN,DAZM, 
we DIST) 

(pes ri = (3N * DISTM) / UNCOV 

Sree soe AC{UDI STM,SE M) 


SLAVE 
Pbieeciven| (A;r ,REOST ,RPOSTL,RSLAV1, RSLALI,FAZ1,BAZ1, 
ep sS TI) 
oesrt = (RN** DIST!) / UNCOV 
ean SSCFAC (UDIST1,SF 1) 

SLAVE2 
meee) INVER? (A,F ,RPOSE,RPOSTL,RSLAV2, RSLAL2, FAZ2Z, BAZ2, 
leo TST 2} 


Coes re 25= (RN * DIST2Z) / UNCOV 
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CY rer Cr] 


Cra YO) 


AOA aA 


200 


Clipe onc r AC {UDIST2Z,SF 2) 


DETERMINE THE ADDITIONAL SECONDARY CORRECTORS FROM THE 
LORAN=-C CORRECTION TABLE FOR THE WEST COAST CHAIN 9940 


CALL TABLE(RPOST,RPOSTL,ID1,ASF 1) 
CALL TABLE(RPOST, RPOSTL,ID2,ASF2) 


GOmeUTE THE RATES AND COMPARE [TO THE OBSERVED RATES 


eso or te UbDTST MY + SF1—- SEN + ASF + DELAY 
Dire’ = TDM1 - RATE! 


iter St Ze DL St meer SF2 = SFM + ASr2 + DELAY2 
Pier 2 = TDM2 = RATEZ 


COMPUTE THE LANE WIDTH IN METERS BASED ON ZQUATION 
4.20 IN RLECTRONIC SURVEYING AND NAVIGATION - 
LAURILA, PAGE 94. 


BR1 = DABS({FAZ1 - FAZ4) 
BR2 = DABS(FAZ2 - FAZM) 
WIDTH1 = (DIFF1 * UNCOV * 0.5) / DSIN(BR1 * 0.5) 
WIDTH2 = (DIFF2 * UNCOV * 0.5) / DSIN(BK2 * 0.5) 


Abel he POSETEON OF VESSEL, COMPUTED RATE, OBSERVED 
RATE, AND THE DIrF BETWEEN THEM. 


Aeon cucu LLALT, IMI N, SSSGC,1 LON, LLYIN, RLSEC,RATE1, 
Peer, RAT EZ, DIF F2 
eee, oem, 0 oo, 5,05, 1X,F6.3,F12.2,F8.2, 


te 12.2, FS. 2) 


WRITE (6,201) 
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QaeY OY 46} 


550 


600 


XMEAN1T AND XMEAN2 ARE THE MEAN DIFFS BETWEEN THE 
COMPUTED RATE AND THE OBERSERVED. XSLAV1 AND XSLAV2 
ARE THE STORED DIFFS. 


XMEAN1 = XMEAN1 + DIFF1 
XMEAN2 = XMEAN2 + DIFF2 
XMEAN3 = XMEAN3 + WIDTH? 
XMEANG = XMEANY + WIDTH2 
XSLAV1(I) = DIFF1 
¥SLAV2(1) = DIFF2 
CONTINUE 


COMPUTE THE MEAN AND STANDARD DEVIATION 


XMEAN1 = XMEANT / IREC 
XMEAN2 ="“XMEAN2 / IREC 
XMEAN3 = XMEAN3 / IREC 
XMEAN4 = XMEAN4 / IREC 


DO 600 I = 1,IREC,1 
VAR1 = VAR1 + ((XSLAY1(I) - XMEAN1) **2) 
VAR2 = VAR2 + ((XSLAY2(I) - XMEAN2) **2) 


VAR1 = VAR1 / (IREC - 1.0) 
VAR2 = VAR2 / (IREC - 1.0) 
SD1 = DSQRT (VAR 1) 
SD2 = DSORT(VAR2) 


WRITE (6,201) 
WRITE (6,210) XMEAN1,SD1,XMEAN3 
WRITE (6,201) 
WRITE(6,211) XMEAN2,SD2,XMEAN4 
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WRITE (6,202) 
210 FORMAT(1X," SLAVE #1, MEAN = ',F10.3, 
1 * STANDARD DEVIATION = ', 
2 F10.3,' DISTANCE IN METERS = ',FP10.3) 
211 FORMAT(1X," SLAVE #2, MEAN = ',F10.3, 
1 * STANDARD DEVIATION = ', 
2 F10.3,' DISTANCE IN METERS = ',F10.3) 


STOP 
See ee eS SSS SS SS SSS SSS STS SSS SS SS SSS SSS SS SSS SSS SSsSsSqeSqsaqzez=e=e= 
SUBROUTINE SECFAC(UTDIST, SF) 
perenne 
C 
c THIS ROUTINE WILL COMPUTE THE SEA SECONDARY FACTOR 
C Tipiste= TONAL DiS TANCE 
C SF = SECONDARY FACTOR 
c 
IMPLICIT REAL*8 (A-H, 0-2Z) 
C 
€ 
€ COEFFICIENTS 
C 
AQ = 129.04323 
A1 = -0.40758 
AZ = 0.00064576813 
G 
Bee. 741282 
B1 = ~0.0711402 
B2 = .000327748 15 
C 
Ma TpisST oGT. 537) GO TO 110 
See ( 80) 7 UTDIST) + B1 + ( B2 * UTDIST) 
GO TO 20 
10 SF = ( AO / UTDIST) + A1 + ( A2 * UTDIST) 
¢ 





Cy OY OP OO OY OY OO OO A AO ae er a a GS 


20 CONTINUE 
RETURN 
END 


ee ee ce es ce es ce ce ee es ee ee es es ee ee ee cece ee ce eee ee ee es ce ee ce es eee es ee ese eee 
er ee ees es es ee ee es es ce es ee ee es ee es ee ee ees es es eee ee eee ee eee eee 


=_—n~pa awe 22 OP eee awe ewe ew ee ee eee ee ee eee Pe eee Pee ee eee eee eee eee eee ee See ee Se ee eee eee eee eee ee eee ee 


SUSKOULENE TABLE SELE CIS THE PROPER ASF CORRECTOR FROM 
imireURAN=C CORRECTION TABLE PUBLISHED BY THE DEFENSE 
MAPPING AGENCY. 
RoAto=— POST EON LATITUDE IN S#CONDS 
RLON = POSITION LONGITUDE IN SECONDS 
LD = DORAN=CechA IN IDENT IPIER 
YW = 1 
= 2 
Y = 3 


Noo DOLE TIONAL SECONDARY «FACTORS 


DHE FOLLOWING TABLES CF ASF CORRECTORS ARE FOR 
Meiers + BAY, CALIFORNEA - S940 -W, -XK, ~Y. 


Pie oeclT REALS (A-H, O-Z) 
PevENohON TASPEW{(3,5) ,TABLEX (3,5) , TABLEY (3,5) 


DATA TADLEW/ 1.6D0,1.6D0, 0.000, 
leo bO, i. 200, 12400, 
fe 200A. 4 OU, OL U, 
1. 4D0,1. 300,1.5D0, 
1. 3D0,0. ODO, 0.0D0/ 


eH WwW NH = 


DATA TABLEX/ -0.9D0,-0.9D0, 0.09D0, 
1 -0.9D0,-0.9D0,-0.8D0, 
-1.0D0,-1.0D0,-1.9D0, 
=-1.0D0,-1.2D0,-1.1D0, 
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OF IY AY aN OF A A 


—& Ww WN 


-1.1D0, 0.0D0, 0.0D0/ 


Pata taSLer7 0. 2D0,0. 200,0.0D0, 
Ga svOyO sO, Os. 3D0!, 
0. 3D0,0. 2D0,0.4D0, 
0.4D0,0. 3D0,0.6D0, 
0. 5D0,0. ODO, 0.0D0/ 
DATA RHOSEC,PI/2.06264806 2471D95, 3.1415926535898D0/ 


TWOPT = PI ¥* 2.0 


CONVER RLAT AND RLON TO SECONDS 


HLAT = RLAT * RHOSEC 
HLON = TWOPT - RLON 
HLON = HLON * RHOSEC 


STARTING LAT AND LONG FOR ScARCH 
LAT = 37700/00.0 LONG = 122/05/00.0 


Pete giiene LHe ASF COR RECTOR FOR THE LORAN-C COMBINATION 


LATITUDE 
ScAr = 133200.0 
SLON = 439500.0 


RMID = 300.0 
Roper = 150.0 


BO 70 g=1,5,1 

SLAT = SLAT - RMID 
ULAT = SLAT + RDIFP 
VLAT = SLAT - RDIF? 
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10 [a(t . OLA mh. LAT .SE. VLAT) GO TO 15 
15 CONTINUE 


LONGITUDE 


DOusOrt—= 1,35, 1 

SLON = SLON - RMID 

ULON = SLON + RDIFF 

VLON = SLON - RDIFF 

30 Pai hoe ULON SANDS hLON .sT. VLON) GO TO 35 
35 CONTINUE 


DETERMINE ASF CORRECTOR 


Powe 0 20 

TP(ID -EQ. 1) ASP = TABLEW (I,J) 
ier DeeHO. 2) ASF = TABLEX (I,J) 
mearp .EQ. 3) ASF = TABLEY (I,J) 
RETURN 


SUBROUTINE INVER? (a,7 INV, GLAT?,GLON1,GLAT2,GLON2,FAZ, 
1 BAZ,S) 


See Appendix B for subrcutine INVERT. 
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DATA SET EXAMPLE ~ PROGRAM LORTAB 


See SePOSLEETION eae LORAN RATES it ee ae 
LATITUDE LONGITUDE 9940-Y in 994 
(== 5} Cp tS) (psec) (usec) Gusec} (nsec) 


36 43 45.800 121 55 27.160 42789.34 -0.29 16294.04 0.34 
Bom+4 3.400 121 55 32.340 42791.13 -0.18 16293.46 0.30 
Boer 246 1900 127 55 37.390 42793.08 -0.18 16292.73 0.41 
Beet 7.490 121 55 46.950 42795.13 -0.38 16292.03 0.34 
meme oss 260 121 55 57.710 42796.93 -0.31 16291.43 0.13 





PROGRAM ASFSEL 


PROGRAN DETERMINES OBSERVED ASF CORRECTORS BY SCANNING 
DATA AT 1 DEGREE LATITUDE AND LONGITUDE INTERVALS. THE 
Poco CrOnS AKE DETERMINED BY SUBTRACTING THE 
CALCULATED TD USING THE SEAWATER SECONDARY FACTOR FROM 
ROE SOBSERVED- TD. RATES . 


CY TOY OO I ae or Oe 


IMELICIT REAL*8 (A-H, 0-2) 
DIMENSION AXIS (13) ,RF (13) 
DIMENSION ASFCR1(16,26) ,ASFCR2(16,26), 
1 INO1(16,26) , IN02(16, 26) 
DATA RHOSEC, PI, UNCOV, RN/2.062648062471D05, 
1 3.1415926535898D0,29 9. 79 2458D0,1.000338D0/ 
DATA XMEAN1,XMEAN2,VAR1,VAR2/0.00D0,0.00D0, 
1 0.00D0,00.0D0/ 
DATA XMEAN3,XMEAN4/0. 00D0,0.00D0/ 
é 
CCC CC C™ RR RK RK KK KKK KKK KKK KKK KK 
DATA AXIS /6.3782064D0 6,6. 378388  D06,6.377397155D06, 
1 6.37816 D06,6. 378 16D 06,6. 378249145D06,6.378165D06, 
2 6.378166D06,6. 378165 D06,6.378145D06,6.3775634D06, 
3 6.378245D06,6. 378135 0D06/ 
DATA RF/6.3565838 D06,2.97 D02,2.9915281285D02, 
1 2.9825D02,2. 9824716 7427D 02,2.93465D02,2.9825D02, 
2 2.983 DO2,2.983D02,2.9825 D02, 6.3562569D06, 
3 2.983 D02,2.9826D02/ 
C#¥*+***ELLIPSOID OPTION NUMBER 
é 1. CLARKZ 1866 8. MERCURY 


le 





Or Gr Orr OO) 


2. INTERNATIONAL (HAY FORD) 2 MARSHAEL ISLAND 


Sa baoonL 194 4 10. NAVY 8D 

4G“. AND (AUSTRALIAN) Iisa ae RY 

Deer O7 REPERENCE 12. KRASSCWSKI 1940 
6. CLARKE 1880 MOD lceeewiosm yo 72 

Ya Ue. 


CCC CCC 4 KK A eK KK ee KK KK KK KK 


c 
C 
c 


100 


Gi. 1 -)- er 7A) 


202 


201 


105 


CC1-2 = ELLIPSOID NUMBER (K) 
CC3-5 = NUMBER OF POINTS ALONG COAST 


READ (5,100) K,IREC,ID1,I1D2 


FORMAT (1X,12,14,213) 


TWOPI=2.*PI 

A=AXIS(K) 

F=1. /RP (K) 

IF (F.LT.3.D-3) F=(A-1./F) /A 


(IREC) 


READ MASTER AND SLAVE STATIONS POSITIONS 
PieeoinateneCORD TS THE NUMBER OF MASTER AND SLAVE 


STATIONS FOR THE PARTICULAR CHAIN. 


meer = (6,202) 
FORMAT (1H1) 
WRITE (6,201) 
FORMAT(1H ) 


READ (5,105) ILATM,IMINM,RSECM,ILONM, ILMINM, RLSECM 
WRITE (6,105) ILATM,IMINM, RSECM,ILONM, ILMINU, RLSECM 
meme (1, 0s, 1X,l2, 1X , PS. 2,1X,14,1K,I2, 1X, FS.2) 
RMAST = ((IABS(ILATM) * 60 + IMINM) * 60 + RSECM) / 


1 RHOSEC 
IF(ILATM .LT. 0) RMAST = -RMAST 
RMASTL = ((IABS (ILONM) * 60 + ILMINM) 
1 RHOSEC 


113 


*60 +RLSECM) / 





PewevONMeeGl-s J) RMWASTL = TWOPI - RMASTL 


Renee tHE FERS@ SLAVE STATION AND CHANGE THE LATITUDE 
AND LONGITUDE INTO RADIANS. 


Ge aa =~ 


READ (5,101) ILATS,IMINS,RSECS,ILONS,ILMINS, RLSECS, 
1 DELAY1 
WRITE (6,101) ILATS,IMINS,RSECS, ILONS, ILMINS, RLSECS, 
1 DELAY1 

meme Homimacy1x,13, 1X,12,1X ,F5-2,1X,04, 1X, 12, 1X,F5.2,F9.2) 
RSLAV1 = ((IABS (ILATS) * 60 + IMINS) * 60 + RSECS) / 


1 RHOSEC 

TE(EWATS -LT. 0) RSLAV1 = -RSLAYV1 

RSLAL1 = ((IABS (ILONS) * 60 + ILMINS) * 60 + RLSECS) / 
1 RHOSEC 


Mm hONo sot. J) RSLAL?Y = TWwOPI - RSLAL' 


c 
C READ THE SECOND SLAVE STATION AND CHANGE THE LATITUDE 
e AND LONGITUDE INTO RADIANS. 
c 
READ (5,101) ILATS, IMINS,RSECS,TLONS, ILMINS, RLSECS, 
1 DELAY2 
WRITE (6,101) ILATS,IMINS,RSECS,TLONS, ILMINS, RLSECS, 
1 DELAY2 
RSLAV2 = ((IABS (ILATS) * 60 + IMINS) * 60 + RSECS) / 
1 RHOSEC 
IF (ILATS .-LT. 0) RSLAV2 = -RSLAV2 
RSLAL2 = ((IABS (ILONS) * 60 + ILMINS) * 60 + RLSECS) / 
1 RHOSEC 


Pee cONnoeeGi. J) XSBALZ = TWOPI ~- RSLAL2 
WRITE (6,201) 


Memo RECORDS POSE TIONS AND OBSERVED RATES, THEN 
SC@yeute THE TOTAL DISTANCE. 
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Domo) I= 1,21 REC ,1 
READS, 104) TELAT, IMEIN ,RSEC,ILON,ILNIN,RLSEC, 
1 RATE, RATE2 

104 EOneQGit, 65,053,613, 15,13,5F653,3%,2F10. 2) 


C CONVERT LAT AND LONG TO RADIANS. 
e 
RPOST = ((IABS(ILAT) * 60 + IMIN) * 60 + RSEC) / 
1 RHOSEC 
IF(ILAT .LT. 0) RPOST = -RPOST 
RPOSTL = ((IABS(ILON) * 60 + ILMIN) * 60 + RLSEC) / 
1 RHOSEC 


IF (ILON .GT. 0) RPOSTL = TWOPI - RPOSTL 


e 
G COMPUTE DISTANCES AND AZIMUTHS FROM THE OBSERVED POINT 
C : 
Cc MASTER 
Ee 
CALL INVER1(A,F,RPOST ,RPOSTL, RMAST,RMASTL, FAZM, BAZM, 
1 DISTM) 
UDISTM = (RN * DISTM) / UNCOV 
CALL SECFAC (UDISTM,SFM) 
C 
SLAVE1 
é 
eC TUENVE RGA, hp RPOST , REOSTL,RSLAVG,RSLAL1,FAZ1,BAZ1, 
1 DIST1) 
UDIST1 = (RN * DIST1) / UNCOV 
Sab SECFAC (UDI ST1,SF 1) 
Cc 
e SLAVE2 


(oO 
a 
CQ 
NO 
Sey 


CavrernvEmh{A,F ,RPOST ,RPOSTL,RSLAV2Z, RSLALZ, FAZZ, 


lat 





1 DIST2) 
UDIST2 = (RN * DIST2) / UNCOYV 
CALL SECFAC (UDIST2,SF 2) 


GS 
ce DETERMINE WHICH LAT AND LONG THE ASF CORRECTOR 
c iS ASsSsi1GNED TO. 
C 
SAUL PASSE GN(REOST APO STLI INI, IN2) 
COMPUTES THE RATES AND COMPARE TO THE OBSERVED RATES 
€ 
Mae — UDtor 1 = UDISTM + SFI - SFM + DELAY 
DIFF 1 = TDM1 - RATE 
¢ 
Poe UNGSotT2 —- UDISTN + S¥F2 - SFX + DELAY 2 
Dees 2e- IDM2 — RATEZ 
eC 
C SUM THE DIFFERENCES TO THE MATRIX AND COUNT THE NUMBER 
e Wiese sCOmnec TORS FOR EACH BLOCK TO LATTER DETERMINE 
c TH MEAN. 
C 


ASFCR1(JN1,JN2) = ASFCR1(3N1,5N2) + DIFF 
ASFCR2(JN1,JN2) = ASPCR2(JN1,5N2) + DIFF2 
INO1(JN1,JN2) = INO1(JN1,JN2) + 1 

550 INO2(JN1,3N2) = INO2(JN1,3N2) + 1 


c Peek Ne THE MEAN AS® CORRECTOR FOR EACH LAT AND LONG 


DO 650 I=1, 16,1 

WRITE (6,201) 

DO 600 J=1,26,1 

IF (INO1(I,J) -.EQ. 0) GO TO 580 

MoGGRIGE,d) = ASECR1(1,0) / IND1(i,J) 
pieee LF (INO2 (I,J) .EQ. 0) GO ro 590 
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Momence (fj) = ASFCR2 (2,3) / IENO2{T,J) 
See ES (ASPORI (1,3) -£0. 9.0) ASFCR1 (I,J) 
eer IF (ASFCR2 (I,J) .£O. 0.0) ASFCR2 (I,J) 
650 CONTINUE 


Je99 
3599 


c 
e WRITE THE CORRECTORS IN MATRIX FORMAT 
c 
C 
IF(ID1 .EQ. 1) WRITE(6, 1000) 
IF(ID1 .£Q. 2) WRITE(6,1001) 
IF(ID1 .EQ. 3) WRITE(6, 1002) 
WRITE (6,201) 
E 


WRITE (6,300) ((ASFCR1 (I,J) ,I=1,16) ,J=1, 26) 
BO 700 L=1,5, 1 
WRITE (6,201) 

201 FORMAT(1H ) 

700 CONTINUE 


IF (ID2 .EQ. 1) WRITE(6, 1000) 
IF(ID2 .EQ. 2) WRITE(6,1001) 
TF(ID2 .EQ. 3) WRITE(6, 1002) 
WRITE (6,201) 
WRITE (6,300) ((ASFCR2 (I,J) ,I1=1,16) ,J=1, 26) 
300 FORMAT(16F6.2) 
1000 FORMAT(' TABLE FOR 9940-W ') 
1001 FORMAT(' TABLE FOR 9940-x ‘) 
1002 FORMAT(' TABLE FOR 9940-¥ ') 


STOP 
END 

C2225 22225 S25 22 eee SH He See se = ee ee ee ee ee = SF See ee eee = = == == ===== 
SUBROUTINE SECFAC(UTDIST,SF) 

Se OOS Oe OO OO Bee BOs 2B 2e eee O28 Oe a Bees es Bee ee eo eee eee @ ee woe wee oe 
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MYA A YN 


OA Aa 


10 


THIS ROUTINE WILL COMPUTE THE SEA SECONDARY FACTOR 


(pts t —~ TOTAL DES TANCE 
SF = SECONDARY FACTOR 


IMPLICIT REAL*8 (A-H, 0-2) 


SOcErICLZENTS 

ROm= 9 129504323 

A1 = -0.40758 

A2 = 0.00064576813 
BO = 2.741282 

eae —-0.011402 

B2 = .000327748 15 


IF(UTDIST .GT. 537) GO TO 10 

SF = ( BO / UTDIST) + B1 + ( B2 * UTDIST) 
50) TO 20 

SF = ( AO / UTDIST) + A1 + ( A2 * UTDIST) 


CONTINUE 
RETURN 
END 


QlFt@? GCYrvG) “Gira 


SUBROUTINE ASSIGN(RLAT,RLON,I,J) 


_e_aap ae SP RP awe ae aa eee eee eee eee eee eee eee ee Se ee manpnpanpPanDp =e eaPRaaPReeapP aH we =weameaweaeae Se Se Se SS ee SP eae ase as = =a» 


PUSROUPENE ASSIGN SELECTS THE COLUMN AND ROW FOR 
PateoalelLUDS AND LONGITUDE OF THE RECORD. 

RLAT POSeeetON LATEIUDE IN SECONDS 

RLON BOs wONecONGITUDE IN SECONDS 

ie COLUMN 
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OQ 1 OC) 


10 
15 


J = ROW 


PPE LeLT REARS (A-H, 0O-Z) 
DATA RHOSEC,PI/2.062648062471D05, 3.1415926535898D0/ 


PWOPT = PTs" 2. 0 
CONVER RLAT AND RLON TO SECONDS 


HLAT = RLAT * RHOSEC 
PEON = TWOPT] = KRLON 
HLON = HLON * RHOSEC 


STARTING LAT AND LONG FOR SEARCH 
LAT = 37/05/00. 0 LONG = 122/05/00.0 


133200.0 
439500.0 


SLAT 
SLON 


LATITUDS 


J = 0 

DO 10 IC = 1,26,1 

Some = SLAY - 60.0 

RULAT = SLAT + 30.0 

RELAT = SLAT - 30.9 

MmtroAdie.s Gl. RULAT sAND. HLAT .GE. RLLAT) GO TO 15 
J = 1c 


LONGITUDE 

I = 0 

Be=20 IC = 1,16,1 
SoeN = SLON - 60.0 


RLLON = SLON + 30.0 


iS 





RRLON = SLON - 30.0 
20 PcOveewnt. RLEON -AND. HLON .GE. RRLON) GO TO 25 
25 r= se 

RETURN 

END 


SUBROUTINE INVER1(A, FINV, GLAT1,GLON1,GLAT2,GLON2,FAZ, 
es A Zig) 


=_ ee ae ae en eee eee SEP ee eee lee ee ee ee eee See eee ee eee eee SP ee eee ee eee ee eee ee ee ee See eee eer ee ee ee eee See ee er ee eee ee eee SP 


See Appendix B for subreutine INVER]. 
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Devas ott EXAMPLE - PROGRAM ASYFSEL 


eyeLE FOR 9940-Y 


1229 / QO! 

Bot 

=-0.47 

=Ve03 “0.44 

Le * 0. —Oe0 | -0.46 -0.40 

en ue. So — 0.5) 0.45 

SOe0S =Oeo7 ——0 549 

=-Q.04 =0.50 

=O 1072 
4or 
ED 
369 
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et oc or 





TABLE FOR 9940-Y 


a0" 


us 


4or 


a! 


Hea 7 55." 
-0.45 -0.47 
Oe) 1 —O0. 50 
mUeoor— 0. 49 
-0.54 -Q.51 
e+ =O. 5/ 
mmo oe = 0s o> 
mu, o? —0.6% 

-0.64 


Seedy 1 ae, 


U0 = 0.40 
ries 
-0.49 
Ue D0 
=U o 
=i. 


Vane 50" 
Omar 
~-90.43 -0.44 -0.47 -90.44 
=Oeoe 60. 46 


Ze 
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